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We investigate the low-energy properties, especially the low-energy excitation structures, of A'^-leg 
integer-spin ladders and tubes with an antiferromagnetic (AF) intrachain coupling. In the odd-leg 
tubes, the AF rung coupling causes the frustration. To treat all ladders and tubes systematically, 
we apply Senechal's method [Phys. Rev. B 52, 15319 (1995)], based on the nonlinear sigma model, 
together with a saddle-point approximation. This strategy is valid in the weak interchain (rung) 
coupling regime. We show that all frustrated tubes possess six-fold degenerate spin-1 magnon bands, 
as the lowest excitations, while other ladders and tubes have a standard triply degenerate bands. 
We also consider effects of four kinds of Zeeman terms: uniform, staggered only along the rung, 
only along the chain, or both directions. The above prediction of the no-field case implies that 
a sufficiently strong uniform field yields a two-component Tomonaga-Luttinger liquid (TLL) due 
to the condensation of doubly degenerate lowest magnons in frustrated tubes. In contrast, the 
field induces a standard one-component TLL in all other systems. This is supported by symmetry 
and bosonization arguments based on the Ginzburg-Landau theory. The bosonization also suggests 
that the two-component TLL vanishes and a one-component TLL appears, when the uniform field 
becomes larger for the second lowest magnon bands to touch the zero-energy line. This transition 
could be observed as a cusp singularity in the magnetization process. When the field is staggered 
only along the rung direction, it is implied that the lowest doubly-degenerate bands fall down with 
the field increasing in all systems. For final two cases where the fields are staggered along the 
chain, it is showed that at least in the weak rung-coupling region, the lowest-excitation gap grows 
with the field increasing, and no critical phenomena occurs. Furthermore, for the ladders of the 
final two cases, we predict that the inhomogeneous magnetization along the rung occurs, and the 
frustration between the field and the rung coupling can induce the magnetization pointing to the 
opposite direction to the field. All the analyses suggest that the emergence of the doubly degenerate 
transverse magnons and the single longitudinal one is universal for the one-dimensional AF spin 
systems with a weak staggered field. 

PACS numbers: 75.10. Jm,75.40.Cx,75.50.Ee 



I. INTRODUCTION 

Low dimensional quantum spin systems have provided 
much interest for a long time. In particular, the un- 
derstanding of one-dimcnsional (ID) spin-^ systems has 
shown a significant progress. Recently, quasi-lD systems, 
such as ladders and tubes, have been among the central 
issues. Here, spin tubes means cylinder-type spin sys- 
tems, i.e., spin ladders with the periodic boundary con- 
dition along the rung (interchain) direction. 

In the spin ladders with an antiferromagnetic (AF) 
intrachain coupling, one of the most dramatic proper- 
ties is the following "even-odd" nature, which is an ex- 
tension of the Haldane conjectur o^i^i'^ for the single AF 
spin chain. For odd- leg and half- integer-spin cases, there 
exist massless excitations above the ground state (GS), 
the spin correlation functions decay algebraically, and 
the low-energy physics is described by a one-component 
Tomonaga-Luttinger liquid (TLL),'* which is equal to a 
conformal field theory (CFT)^ with the central charge 
c = 1. Meanwhile, for other (even- leg or integer-spin) 
cases, the system is gapful, and the decay of spin corre- 
lations is an exponential type. This even-odd property 
has been established by both numericalSi^iSiS and analyt- 



icaU&^^i^ works. Moreover several experiment3i2iiaiSii& 

also support it. 

Both theoretical and experimental studies of spin tubes 
are not as active as those of spin ladders. As far as 
we know, there are only two spin-tube-like materialsi^iiS 
even now. However, odd-leg tubes with an AF rung cou- 
plingi9-20-2i-22,23.24,25 i^^yg attracted considerable inter- 
est at least theoretically, because such tubes possess the 
frustration along the rung. It is known that at least for 
the strong rung-coupling regime, odd-leg AF-rung spin- 
i tubes take doubly degenerate and gapful GSs with the 
one-site translational symmetry along the chain breaking. 
Namely, the rung frustration induces the break down of 
the even-odd prediction. 

As powerful theoretical tools to treat these ladders 
and tubes, there are nonlinear sigma model (NLSM) ap- 
proachesiiiii2*2SiSL22i2a A standard NLSM technique, 
which has an ability to derive the above even-odd na- 
ture, assumes the development of a sufficient short-range 
order to all spatial directions. Thus, it is not applicable 
for frustrated odd-leg tubes. However, if we first map 
a single AF spin chain to a NLSM, and next take into 
account the rung coupling perturbatively, we can deal 
with frustrated tubes as well as other non-frustrated sys- 
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tems within the NLSM framework. In this paper, follow- 
ing this idea, we revisit and investigate the low-energy 
physics of spin ladders and tubes systematically, in the 
weak rung-coupling regime. Note that the perturbative 
treatment of the rung coupling was already proposed by 
Senechal, who applied it to 2-leg ladders. Therefore, 
our method discussed below will be regarded as a nat- 
ural extension of his work. As well known, the NLSM 
method for half-integer-spin chains bears a topological 
term (Berry phase). Because (as Senechal mentioned 
in Ref. Isol) it is difficult to treat such a term and the 
rung coupling concurrently, we concentrate on integer- 
spin cases only in this paper. 

Our target is the following Hamiltonian for A^-leg spin 
systems: 
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N 



n 



Sl] ■ Sij+i + J±^^ Sij ■ Si 



(1) 
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where Si,.j is the integer-spin-5' operator on site (Z,j), 
J(> 0) is the intrachain coupling, and J_l is the rung 
one. In the rung-coupling term, ladders take N = N — 1, 
while N = N and Sn+i,j = Sij in tubes (TV > 3). 

We further study external-field effects for the 
model In this paper, we consider following four kinds 
of the Zeeman terms: 

^[0,0] = ^H^^i.j^ (2a) 

1,3 

= -HY,i-^y+'Sl^, (2b) 

1,3 

H[.,o] = -HY,i~^ySf^,, (2c) 

1-3 

= -i/5](-l)'+-''+^5,:,, (2d) 

1,3 

where H{> 0) is the strength of the external field. 

The first term 7i[o.o] is a standard uniform-field Zee- 
man term. External fields of other terms have 
an alternation. We call the fields in H2a() - (|2d|) 
as a [0,0] (uniform), a [0,7r] (staggered along the 
rung), a [tt, 0] (staggered along the chain) and a 
[n,TT] (staggered along both directions) fields, respec- 
tively. Staggered magnetic fields have been investigated 
recently 3ii32,33,34,35,36,37,38,39,40,4i^42,43,44 Actually such 
fields are present in real magnetS)36i37i45,46,47,48,49,5Q,5i 
and their origins have been explainedi^SiiS One will see 
that these four terms are congenial to the NLSM method. 

The organization of this paper is as follows. First, we 
review the NLSM approach for single AF integer-spin 
chains in Sec. |nj It provides an underlying effective the- 
ory for spin ladders and tubes. Section UTTI presents our 
main results, in which we treat spin ladders and tubes 
in quite detail. Sections ITTra UTrei UTTO and UTTdI are 
devoted to investigate the no-field, uniform([0, 0])-field, 
[0, 7r]-field, and [tt, 0] or [tt, 7r]-field cases, respectively. 



One will see several new even-odd natures in spin lad- 
der and tube systems. Particularly, in the no-field and 
[0, 0]-field cases, we find qualitative differences between 
the low-energy excitation in even-leg tubes and that in 
odd-leg (frustrated) ones. Because these results in the 
two cases are supported by symmetry arguments, we be- 
lieve that they are not merely approximate results, and 
true. In Sec. IIVI we summarize all the results and touch 
some related topics. We write down the properties of 
some simple matrices in Appendix ^ Moreover, Ap- 
pendix Ogives a review of Green's functional treatment 
of the staggered field along the chain. These Appendices 
are useful for the calculations in Sec. IIIII 



II. REVIEW OF SINGLE-CHAIN CASES 

We present a review of the NLSM approach for single 
integer-spin AF chains and the saddle-point approxima- 
tion (SPA) in this section, which mainly follows Refs.l39l 
and El. This will be the basis of Sec. |llll Readers 
who are familiar with the approach can skip this section, 
especially Sec. Ill Al 



A. No-field case 

This subsection discusses the integer-spin-S" isotropic 
Heisenberg AF chain without external fields. 
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There exist two celebrated ways to obtain the low-energy 
effective theory for the chain Q, a NLSM: the operator 
formalismii*2& and the path-integral oneiSLSLSS We use 
the latter here. In the latter formalism, the Euclidean 
action of the chain Q is given by = dTTi[(i{T)], 
where r is the imaginary time, and /3 = l/(fcBT) is the in- 
verse of temperature. The quantity H[r2(r)] is the "clas- 
sical" Hamiltonian in which the spin operator Sj/S is 
exchanged into a three-component unit vector f2j (r) . 

Following Haldane's ideaji we take the spatially con- 
tinuum limit Qj (r) — > Q(xj , r) [xj = jxa: a is the lattice 
constant], and decompose Q{x,t) into the uniform fluc- 
tuation 1{x,t) and the AF fluctuation n{x,T) as follows: 



(i{xj) w {~iyn{xj)\Jl - a'^l{xjy + al{xj), 



(4) 



where two new constraints fp — 1 and n ■ I — are 
imposed to maintain the original constraint il^ = 1 up 
to 0{a). The approximation |0J is called the Haldane 
mapping, and it assumes that there exist the low-energy 
and slowly-moving modes around both the wave num- 
bers fc = and k = ir/a. It is hence expected that the 
more the GS approaches a Neel ordered one (i.e., the 
classical limit S — > oo), the more the mapping is reliable. 
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Through a few procedures [(i) substituting Eq. to Ab, 
(ii) a gradient expansion for A^, and (iii) integrating out 
the uniform part I or replacing / with its classical solution 
Ici = 45 ja ("^ X which is defined by 5A-e,/51 = 0], we 

can finally obtain the effective model for the AF fluctu- 
ation 71, 



Z 



j T>fiD\ cxp(— 5e[7i, A]), 



dx. £e(^(x), A(x)), 



(5a) 
(5b) 



2g 



dl + cdl 



i\{rf - 1), (5c) 



where 5e, and £e are the partition function, the Eu- 
clidean action, and the Lagrangian density, respectively. 
The symbol x means (x, t), g = 2/ S is the bare coupling 
constant, c — 2SJa is the bare spin- wave velocity, and 
A(x) is the auxiliary field for the constraint = 1. The 
model (O is nothing but an 0(3) NLSM. In this frame- 
work, the spin operator is approximated as 



Sj w (-iySn+—{nxn). 



{n = drfi) (6) 



Using this model, let us consider the low-lying band 
structures in integer-spin chains. It has been known well 
that the (l-M)D 0(3) NLSM is integrable:'^^ the system 
is gapful and the first excitation bands consist of the 0(3) 
triplet particles. However, the integrability method can 
not be extended to the case of ladders and tubes 
In this paper, we utilize a saddle-point approximation 
(SPA) instead to reveal the low-energy properties quali- 
tatively, in the simplest manner. (As one will see later, 
the SPA is available in ladders and tubes.) Integrating 
out the field 7?(x) in Z, we obtain the effective action 
including only the field A(x), 5e[A], which is defined by 
Z = / I?Ae"'^'=[^l. The SPA in the present work is given 
by replacing A(x) with Agp (a constant independent of x 
and r) which is the solution of the saddle-point equation 
(SPE) 95e[Asp]/9Asp = 0. To obtain the exphcit form 
of S'e[Asp], we introduce the Fourier transformation for 
n"(x) [a — X, y, z] as 



n"(x) 



1 



dx \ dr e 





X — lUJnT-OL 



-''=^+■"^"^71" (x), (7a) 
n"(a;„,fc), (7b) 



where L = Ma is the system length (M is the total site 
number), k = i'Km/L is the wave number, and LOn = 
27m//3 is the Matsubara bosonic frequency (rn^n G Z). 
Hereafter, we will often use a new symbol k = (a;„,fc). 
Because n(x) is real, Ti"(k)* = n"(— k). Performing the 
Gaussian integral of the field fi"(k) in Z, we obtain 



^eIAsp] 



2 ^ 

k 



In 
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{uji + c'k') - iX 



-iL(3\sp. 



(8) 



Therefore, the SPE is evaluated as 



3gc r dk , //? ,,,\ 



(9) 



where e{k) = c^/k^ + ^^, = — i2gAsp/c and A is the 
ultraviolet cut off. We performed the sum of ti>„ using the 
standard prescription with the residue theorem,^'* and 
then took the continuous limit ~* / if - limit 
T ^ 0, Eq. ® is reduced to 



3gc 
2tt 



dk 



27r 



-1- (10) 



Equations ^ or H1U|) fix Agp and ^. They also suggest 
that ^ is real and Asp is purely imaginary. To complete 
the SPA, we must determine the unknown parameter 
A. In other words, the present SPA provides a one- 
parameter-fitting theory. Of course, other quantities such 
as g and c can be adopted as the fitting parameter. How- 
ever, we will always take A as it throughout this paper. 

The SPA transforms the constraint to a mass term 
for the bosons n". Therefore, after the SPA, the field 
n(x) stands for the triply degenerate massive bosons with 
dispersion e(fc). This is consistent with the exact solution 
of the NLSM. The bosons should be regarded as the spin- 
1 magnon excitations of integer-spin chains Q . The gap 
e(0) — c^^^ hence corresponds to the Haldane gap A. 
From Eq. H1U|I . we obtain 



e(0) = cA/sinh(S'7r/3) , 



(11) 



It is remarkable that the Haldane gap depends on the 
spin magnitude S in an exponential fashion. It corre- 
sponds to the fact that the conventional spin-wave the- 
ory (1/5' expansion) can not explain the Haldane gap. 
The accurate values of A in spin-1, 2 and 3 AF chains 
are found by numerical works Therefore, we can de- 
termine the cut off A from the relation A — e(0), com- 
pleting the SPA. The NLSM plus SPA scheme further 
leads to {S'^S^) « {-ly S^{n'^{xj)n°'{Q)) + ..., {n^) = 1 
and (n^'^'^) = ((• • • ) stands for the expectation value). 
The first result means that ^ is interpreted as the spin 
correlation length. The second is trivial from the SPE, 
= 9iSe/9Asp oc {"n? — 1). The third is also trivial due 
to invariance of the action Q under n —n. 

Table m provides the numerical data [quantum Monte 
Carlo (QMC) simulation, exact diagonalization and 
density-matrix renormalization-group (DMRG) method] 
and the above SPA results for the spin-1, 2 and 3 chains. 
From this, the SPA is expected to work well even in 
the minimum-integer-spin (spin-1) case. The larger the 
spin magnitude S becomes, the more the SPA correlation 
length ^ approaches its correct value. This is consistent 
with the fact that the NLSM method is considered as an 
expansion from the classical limit (S — *■ oo), i.e., a Neel 
state. On the other hand, the effective Brillouin zone (or 
the cut off A) rapidly becomes smaller with increasing S. 
This implies that the SPA is efficient only for extremely 
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TABLE L Haldane gaps, spin-spin correlation lengths and spin-wave velocities in spin-1, 2, and 3 AF chains @. Th e Q MC 
data of Haldane gaps and correlation lengths are quoted from Ref. IbrL The velocity c of numerics is found in Refs. l57l and 





Haldane gap A 




correlation length ^ 


spin-wave velocity c 


Spin 


(QMC data) 


optimized cut off A 


QMC data 


SPA data 


numerics 


bare value 


1 


0.410 X J 


0.0816 X n/a 


6.015 X a 


4.872 X a 


^ 2.5 X Ja 


2 X Ja 


2 


0.0892 X J 


0.02837 X TT/a 


49.49 X a 


44.858 X a 


^ 4.65 X Ja 


4 X Ja 


3 


0.0100 X J 


0.006139 X n/a 


637 X a 


598.80 X a 


? 


6 X Ja 



low temperatures; fcsT ^ A. Of course, one can con- 
tinue more precise analyses of the NLSM Q beyond the 
SPA (for example, using its exact solution, renormaliza- 
tion group, large-iV expansion, the improvement of the 
magnon dispersion, etc) iSSiStSiffi 



B. Uniform-field case 

We consider integer-spin chains ^ with the uniform 
Zeeman term, —Hj^j^j- Recalling that the boson n 
in the NLSM represents the triply degenerate spin-1 
magnons in the no-field case, one can immediately con- 
clude that the uniform field splits the degenerate bands 
into three ones, which have 5^ = 1, 0, and —1, respec- 
tively. In this subsection, we verify that the NLSM and 
the SPA can reproduce this Zeeman splitting. 

Within the NLSM formalism, the uniform field H = 
{0,0, H) couples to the uniform fluctuation Im^ In this 
case, the classical solution of I becomes 



magnon dispersions. Therefore, the field H induces the 



ASJa 



[fi X n) 



1 



45Ja 



(H + iXn), (12) 



where the second term in the right-hand side originates 
from the field H, and A is the auxiliary field for the con- 
straint I ■ n ~ 0. We fix A to iH ■ n which is the solution 
of Ici ■ n = 0. Substituting Id to the low-energy action, 
we obtain the action of n. 



H ■ (n X n) 

4Ja ^ '. 



(13) 



where £e is the same as Eq. (I5c|l . and remaining two 
terms are induced by the uniform field. Since the action 
is also quadratic in the field n, it is possible to integrate 
out it in As a result, the SPE for A is 



27r 



dk X ^ 

eO(fc) ^ 



coth 



(14) 



where e°{k) = c^k"^ ^ S,-^ + H'^/c^ and e^{k) = e°(fc) =F 
H. We use the same cut off A in the no-field case. 
Observing the Fourier-space representation of the ac- 
tion H13|l or calculating the two-point correlation func- 
tions of n, one finds that e°'+'~(fc) are regarded as the 



band splitting e(fc) 
re-expressed as 



.0, + , 



(fc). At T = 0, Eq. (HH) is 



3gc 
2tt 



dk 
^o(fc) 



1. 



(15) 



The comparison between Eqs. (jlOf) and (|15ll shows that 
^{H = 0)-2 i{H)-'^ + TJVc^ and e(fc) = t^{k) are re- 
alized at T = 0. These two relations tell us that the SPA 
reproduces the Zeeman splitting of the spin-1 magnon 
modes at T = 0. Modes e°(fc), e~''(fc) and (~{k) can be 
regarded as = 0, +\ and —1 magnons, respectively. 

Similarly to the preceding subsection, the SPA derives 
(n^) = 1 and (n^'*''^) = 0. However, it also does an 
incorrect result {S^) cx {l^^ ^ 0. It would be because the 
SPA and the Haldane mapping do not take care of the 
spin uniform part / sufficiently, in comparison with the 
staggered part n. 



C. Staggered- field case 

Let us next discuss integer-spin chains lj3Jl with the 
staggered Zeeman term, — H 'Y^ ^{—ly , in which the 

staggered field {—lyH directly couples to the AF fiuctu- 
ation n. The staggered magnetization = (— 1)^ (5*1), 
magnetic susceptibilities and transverse excitations can 
be evaluated by applying the SPA in sufficiently low tem- 
peratures, like Sees. Ill K\ and ITT^ However, it has been 
shown in Ref. [s^ that in addition to the SPA, the Green's 
function method is necessary for a quantitative estima- 
tion of longitudinal excitations. "Transverse" ("longitu- 
dinal" ) means the components of n which are perpendicu- 
lar (parallel) to the staggered field. Here, we provide only 
main results of the Green's function method in Ref. l39l 
Its brief explanation is in Appendix^ which is applied 
in Sec. IIII CI For more details, see Ref. Isol 

The normalized magnetization (H) /S is fixed by 
Eq. HB12|I . We draw it in Fig. ^ which explains that 
as S becomes larger, m^{H)/S grows extremely rapidly. 
The Green's function method (plus SPA) derives two-fold 
degenerate transverse magnons with the dispersion exik) 
and a nondegenerate longitudinal magnon with eL(fc). Of 
course, these two bands return to e(fc) as = 0. The 
transverse gap (lowest excitation energy) At — exiO) 
and the longitudinal one A^ = 6^(0) are determined from 
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FIG. 1: Normalized staggered magnetizations ml / S in spin-1, 
2, and 3 chains with the staggered field at T = 0. 
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FIG. 2: Transverse and longitudinal gaps in spin-1, 2, and 3 
chains with the staggered field. We set J = 1. 



the SPE ||BTT|| and the relation ||B22|| . Figure |3 repre- 
sents two gaps At and A^. It shows that the gap grows 
more rapidly as a function of H for larger S. The rela- 
tion At < At < 2 At always holds within the present 
Green's function nicthodii^ 

In Ref. 39, it has been confirmed that in the spin- 
1 chain, m^, At and At excellently agree with those 
determined from DMRG method within the weak-field 
regime {H < J).''^ Recalling again that the NLSM is 
a semiclassical approach, we can guess that the above 
quantities of the spin-2 and 3 chains, determined by the 
NLSM, are also consistent with their correct values at 
least in the regime H <^ J. 



III. LADDERS AND TUBES 



Based on the NLSM method for the single-chain prob- 
lems in the preceding section, we investigate our targets, 
iV-leg integer-spin ladders and tubes ^ in T = 0. 



A. No-field case 

This subsection considers TV-leg integer-spin ladders 
and tubes without external fields. (Here we remark that 
the 2-leg spin-1 case has already been analyzed by the 
NLSM,^ a non-Abelian bosonization^ and the QMC 
simulation. Following the method in Ref. [s^, we treat 
the rung-coupling term as a perturbation on N decou- 
pled chains, each of which can be mapped to an NLSM. 
Namely, within the NLSM framework, we approximate it 
as follows: 



1=1 J 



Si+ij 



N 

E 

1=1 



dx fii ■ fii+i 

(16) 



where fti is the AF fluctuation field of the l-th chain 
(nN+i = ^i)- This prescription enables us to deal with 
any rung-coupling terms even including frustrations, al- 
though it would be valid only in the weak rung-coupling 
regime; J S> \J±\- However, as a price, we have to take 
into account N constraints: ftf = 1 {I = 1, . . . ,N). (As 
well known, there is only one constraint in the standard 
NLSM method.) Here, as a crude approximation, we re- 
place these constraints with an averaged one. 



N. 



(17) 



We will discuss the validity of Eqs. (|16|l and H17() in the 
final part of this subsection. (As one will see there, we 
can predict that these two approximations are allowed in 
the any-leg weak-rung-coupling systems, at least within 
the qualitative level.) Under the approximations H16() 
and l|17|l . the total action of ladders or tubes is described 
as 



SE[{ni},X] 



(18) 



where the subscript means transpose, Af^ = 
^(n", . . . , n^), A is the auxiliary field for the constraint 
(|17|l . and the N x N matrix A is defined as 



A = 



ai 

02 



25c 

S^J±/{2a), 



ao \ 



■■ 02 

02 ai ) 



iX, 



(19a) 

(19b) 
(19c) 



where ao = (02) for ladders (tubes). It is notewor- 
thy that the action of an AF-rung ladder (an even-leg 
AF-rung tube) can be transformed to that of the ferro- 
magnetic (FM)-rung ladder (the FM-rung tube) through 
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the unitary transformation, rii=even — * ~ni=even- There- 
fore, both AF- and FM-rung ladders (both even-leg AF- 
and FM-rung tubes) with the same leg number have the 
same low-energy excitation structure within the present 
scheme. Indeed, this property has been partially ob- 
served in a QMC study of the 2-leg spin-1 ladder with 
<C On the other hand, we also emphasize that 
there are no unitary transformations connecting an odd- 
leg AF-rung (frustrated) tube and the FM-rung one. 

Because the action lfTH|) is quadratic in the fields n/ 
like the chain cases, we can integrate out and derive 
the SPE for A. After diagonalizing A (see Appendix IXjl 
and performing the Fourier transformations for Hi (x) , the 
action becomes 

5e = ^^(c.2+e.(fc)2X(k)*<(k) 



-iLpNXsp, 



(20) 



where 



e,(fc) 



cA/k^+^r\ (21a) 
2 + 2J_LCOsA:,./(Ja2), (21b) 



< r < 



and fc, = ^ (^) and r = 1,...,7V (-[f] 

[^]) for ladders (tubes: N > 3). The symbol [v] means 
the maximum integer u satisfying u < v. The new field 
TO"(k) is defined by m"(k) = Uri'h'j' (k) , where Uri is 
the unitary matrix diagonalizing A. In the derivation 



of Eq. H20() . we performed the replacement A 



Asp (a 



constant), and assumed that Asp is purely imaginary. In 
tubes, kr means the wave number for the rung direction. 
Following the SPA prescription similar to that in Sec.lTTl 
we easily estimate the SPE, 9S'e[Asp]/9Asp = 0, where 
5e[A] = -ln(/niM,e"^E[{"'>'^l), as follows: 



27r ^ 



dk 
er(fc) 



coth 



(fe.(fc) 



N, (22) 



where we use the same cut off A as that of chains. For 
J_L = 0, Eq. ^ reduces to the SPE for chains ©. 

The representation H2U|) tells us the following several 
low-energy properties of ladders and tubes, (i) One can 
interpret er{k) as a spin-1 magnon dispersion. The band 
splitting e(fc) er{k) is due to the hybridization effect of 
the rung coupling. Each band er{k) is triply degenerate 
correspondingly to S"^ = 1, and —1. (ii) The 3A^-band 
crossing occurs at J± — 0. Obtaining this level-crossing 
picture is an advantage of the NLSM approach>2° Other 
methods, such as the QMC simulation^ and the non- 
Abelian bosonization^SiSliS^ cannot derive it. (iii) In ad- 
dition to the triple degeneracy of S^, tubes (not ladders) 
have an extra degeneracy er{k) — e^r{k). Namely, in 
tubes, only two bands eo(fc) and ejv (fc) are triply degen- 
erate, while all other bands have a sixfold degeneracy. 
Here, note that the original tube is invariant under 
reflection including (or tt rotation with respect to) the 




FIG. 3: Cross section of the 5-leg tube and the reflecting 
operation for the rung direction. 



central axis of the cross section of tubes (Fig. Be- 
cause this operation causes k^ — > — fcr, we can conclude 
that the degeneracy er(fc) — e-r{k) comes from the re- 
flection symmetry and it must be a correct result inde- 
pendent of our approximation scheme, (iv) Noticing the 
contents of (i)-(iii) and the form of er{k), we can show 
the band splitting as in Fig. ^ For any r, er{k) has 
the minimum at fc = 0. Thus we define the gap of each 
band, = £^(0). The true gap Amin of the system 
would be the smallest among A^. For the nonfrustrated 
systems, Amin is always carried by a triply degenerate 
band. Those of AF-rung ladders, FM-rung ladders, FM- 
rung tubes, and even- leg AF-rung tubes are given by A^r, 
Ai, Aq, and Aw, respectively. On the other hand, the 
gap of frustrated tubes is carried by a six-fold degenerate 
band with e w-i (fc) = e w-i (fc). This is a new even-odd 

property in the AF-rung spin tubes. This interesting phe- 
nomena can be intuitively understood as follows. The 
GS in all AF-rung tubes would tend to take a short- 
range Neel order along the rung. Therefore, we guess 
that the lowest excitations on such a GS are around the 
wave number kr = n. Actually, those in nonfrustrated 
tubes always have kr = tt. However, the wave number 
fcr = TT can not be admitted in frustrated tubes. Instead, 
the lowest bands in A^-leg frustrated tubes hence have 
fcr = TT ± , which is closest to tt in all the wave num- 
bers. These bands are just sixfold degenerate. On the 
other hand, the lowest band in the FM-rung case always 
has kr = because of the similar reason; FM-rung tubes 
tend to develop a FM short-range order along the rung. 
All the band with fc^ = are not degenerate, except for 
the degeneracy of the spin-1 triplet. 

We investigate the low-energy excitations more quan- 
titatively by solving the SPE H22|l . As we will see in 
Figs. 151101 all ladders and tubes have a positive-^r solu- 
tion, which means that all magnon bands er{k) are well- 
defined. First, we discuss nonfrustrated systems. Fig- 
ure |31 shows the rung-coupling and spin-magnitude de- 
pendence of gaps in 2-leg ladders. Figure El represents 
the rung-coupling dependence of gaps in iV-leg spin-1 lad- 
ders. Moreover, Fig. [7| provides the lowest gap Amin in 
A^-leg spin-1 ladders and tubes. (Since the band struc- 
ture of a nonfrustrated system with Jj_ is equivalent to 
that of the system with in the present scheme, the 
same result applies to the FM-rung ladders, and to even- 
leg AF-rung tubes.) The former two figures indicate that 
the rung coupling induces rapid rises of bands except for 
the lowest bands. It supports the known result that the 
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(|M'-rung ladder N >2 FM-rung ladder N > ^l] 




3N-fold degenerate 
spin-1 magnon bands 



FM-rung tube N > 3 
even-leg p odd-leg 




.£(N-1)/2=E -(N-1}/2 
£(N-3)/2=e -(N-3)/2 



AF-mng tube N > 3 

odd-leg (frustrated) 



r , . , E|N-3)/2=E -(N-3)/2 
mLfl'^^e|N-l)/2=E -(N-1)/2 



FIG. 4: Band splitting induced by the rung coupling J_i_. 
"Three- or sixfold" means the degree of the band degener- 
acy. 
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0.2 




Spin-1 
Al/ 



1.4r 



Spin-2 
1 .0 ■ Ai / 



Spin-3 



0.05 -i 



A2 



0.1 0.2 0.3 0.02 0.05 0.0002 0.0004 
Jl J± J± 



FIG. 5; Gaps Ai and A2 = Amin in 2-leg AF-rung spin-1, -2, 
or -3 ladders with J — 1. 




0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 
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FIG. 6: All band gaps Ar in A'^-leg AF-rung spin-1 ladders 
with J — 1. 




FIG. 7: Lowest magnon gaps Amin in A'^-leg AF-rung spin- 
1 ladders and FM-rung tubes (non-frustrated systems) with 
J = 1. The gap-reduction speed of the A^-leg tube is larger 
a little than that of the ladder. The 2-leg "tube" means the 
2-leg ladder. 



standard NLSM approach for nonfrustrated systems 
which extracts only the lowest bands, captures the low- 
energy physics in | Jj^ | ^ J. Figure shows that the 
increase of N gradually enlarges splitting-magnon-band 
width. One finds from Figs. Inland [7| that the more N or 
S increase, the larger the decreasing speed of gaps Amin 
becomes around the decoupling point J± — 0. Particu- 
larly, in Fig. |S1 it is remarkable that once one attaches 
an awfully weak rung coupling for a ladder with 5 > 1, 
the gap Aftiin sharply approaches zero. For example, 
when we set (J, Jj_,T) = (1,0.05,0), the SPA predicts 
(Amin,^) = (0.268,0.653) for S" = 1, (0.0143,0.160) 
for 5 = 2, and (0.000459, 0.0458) for = 3 (A is the 
Haldane gap of single chains in Table QJ. These gap re- 
ductions are naturally expected from the consideration 
that the growths of S and N help the GS (a massive 
spin-liquid state) be close to a Neel state, which has a 
massless Nambu-Goldstone mode. 

We next focus on the frustrated spin tubes. Figure |H1 
is the gap structures of AF-rung spin-1 tubes. The 
panel (a) shows that gap reductions of frustrated (odd- 
leg) tubes are considerably slower than those of even-leg 
tubes. It must reflect that the rung frustration obstructs 
the rise of the AF short-range order unlike in the non- 



frustrated systems. While, similarly to nonfrustrated sys- 
tems, the growth of TV prompts the gap reduction in the 
frustrated tubes. It would be a relaxation effect of the 
frustration. We see from the panels (b) that all magnon 
bands, except for lowest one, quickly rise together with 
increasing J± even in frustrated tubes. This may sug- 
gest the possibility to construct an effective theory for 
frustrated tubes, which includes only lowest sixfold de- 
generate bands. 

Finally, we discuss the validity of our strategy in this 
subsection. In order to investigate the integer-spin lad- 
ders and tubes, we took approximations (|16|l and H17|l . 
First, we consider the validity of Eq. H17|l . If we adopt 
the original constraint nf = 1 {I = 1, . . . , N) instead of 
the averaged one (|17|l . the action corresponding to the 
model JQ) is given by 

r ^ 1 

S'^[{ni}, {Xi}] = y 5] [n, ( - — {dl + c'dl) 



iXAfti + iXi 



E 

1=1 



ni ■ ni+i 



}, (23) 



where A; is the auxiliary field for the constraint nf = 1. 
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(a) AF-run^ tub es (b) frustrated tubes 
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Ao/ 
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0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 

Jl J± J± 



FIG. 8: Lowest magnon gaps Amin and band gaps A,, in 
A^-leg AF-rung spin-1 tubes with J — 1. 



At least for the 2-leg ladder system, the SPEs for the 



original constraint nf 



1 turn out to be identical 



with that for the averaged one.^ The deviation between 
original and averaged constraints could appear in 3-leg 
or higher-leg systems within the SPA. Actually, one can 
easily check that the SPEs for the original and averaged 
constraints provide different solutions in the 3-leg lad- 
der. Therefore, averaging H17|l is expected to be invalid 
for large-iV systems. Here, notice that the action (|23f) is 
invariant under the transformation ni — s- fti^Q and A; — s- 
Xi+Q [Q e Z] (n; UN+i-i and A; Xn+i-i) in the 
tube (ladder) systems. These transformations of course 
correspond to the translational operation along the rung 
in tubes {Sij Si+qj) and reflection around the cen- 
tral axis along the chain in ladders {Sij Sn+i-ij), 
respectively. Let us discuss the form of the SPEs, us- 
ing these symmetries. Each original SPE is written as 
dSi,[{Xm}]/dXi oc {nf - 1) = 0, where S'[{Xm}] = 
-ln(/nm2?"^me"'^E[{«™},{^™}]). The expectation value 
{nf) is a function of {A™}; {nf) = /(Ai,-- - ,Ajv)- In 
the tube systems, these facts and the above symme- 
try of the action 123|) lead to the identity 1 = {nf) = 
/(Ai, • • • , Aat) = {nf_^_^) = f{X2,- ■ ■ ,Xn,Xi) = 
Therefore, we find that dS'^[{Xm}]/dXi\x^=...=Xj^ is in- 
dependent of all the chain indices I in tubes. On the 
other hand, in tubes, the SPE for the averaged constraint 
concerns those for the original ones as follows: 



= 



dSEjX] 
dX 



N 



dS^[{X,r,}] 



1=1 



N 



dXi 

dS'j,[{Xrr.}] 



9Ai 



(24) 



where the final equal sign is thanks to the above property 
of the function /. Because we have already known that 
(9S'e[A]/9A = has a physically suitable solution Asp, 
Eq. (|24ll indicates that the original SPEs in tubes can 
take the same solution, Ai = • • • = Aat = Agp. Thus, the 
averaging of the constraints, Eq. H17|) . should be valid 
on the symmetric solutions of the original SPEs, Ai = 



• • • = Aat , although other possible solutions would not 
be covered. Meanwhile, in ladder systems, the similar 
argument can not lead to the validity of Eq. (|17|1 . If the 
effective theory (|23f) for ladders possesses the refiection 
symmetry {Sij Sn+i-ij), the original SPEs should 
have a solution with A/ = X^+i-i- This solution does 
not contradict the symmetric one. Therefore, we expect 
that the averaging (|17|l is admitted even in ladders. 

Subsequently, we discuss the perturbative treatment of 
the rung coupling Hlt)|) . As mentioned before, the approx- 
imation would be justified only in the weak rung- 
coupling regime \J±\ <^ J. As we see from Figs. 
the present scheme always predicts that the gap Amin 
monotonically decreases with |Jj^| increasing in all sys- 
tems. However, it has been known, from the QMC sim- 
ulatioujS that in the 2-leg spin-1 AF-rung ladder, the 
gap reaches its minimum value at a finite Jl, and then 
grows and approaches the gap for the 2-spin problem in 
the single rung with increasing J± . (The simulation also 
shows that the gap monotonically decreases for FM-rung 
side.) In addition, the standard NLSM method for non- 
frustrated systems fii*i2iffi which is reliable for the case 
I Jl I ^ J, shows that the gap is a monotonically increas- 
ing function of J± in all the AF-rung cases. These gap 
growths cannot be explained in our weak rung-coupling 
framework. (Inversely, it also means that the standard 
NLSM approach breaks down in the weak rung-coupling 
regime.) Figure |51 displays gaps Anun of spin-1 ladders 
and tubes with a few leg numbers, which are determined 
from the SPA and the QMC method (note that the QMC 
method cannot be applicable for the frustrated tubes doe 
to the negative-sign problem). As expected, the SPA gap 
is semiquantitatively identical to the QMC one within the 
sufficiently weak rung-coupling regime |J_l| ^ 0.05 x J, 
outside which the deviation between SPA and QMC gaps 
becomes clear. The good agreement between the SPA 
and QMC methods is observed in 2, 3, and 4-leg sys- 
tems. It encourages and allows us to apply the present 
SPA scheme to large- TV systems. The SPA gap is always 
larger than the QMC one in the region | J_l| > 0.05 x J. 
It might be because the SPA does not take into account 
the quantum fluctuation effects enough. In Fig. ^| we 
draw the spin-spin correlation lengths of small-A^ sys- 
tems, which are determined by the SPA and the QMC 
methods. There, we define the SPA correlation length as 
^corr = cA~?jj. Since our scheme optimizes the gaps and 
not the correlation lengths, the deviation of the SPA and 
QMC data already emerges at the starting point Jj^ = 0. 
However, like the gap, the behavior of the SPA correla- 
tion lengths is similar to that of the QMC ones within 
|J_l| < 0.05 X J. 

We believe that our results are qualitatively correct in 
all ladders and tubes with a weak rung coupling, and 
in particular, the predicted band structures in Fig. ^ 
are true. The band degeneracy is strongly protected by 
the symmetry argument. Like the final statements in 
Sec. Ill AI one can imagine several modifications of ap- 
proximations Hlt)|) and l|17|l . 
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4-leg ladder 
and tube 



0.05 , 0.1 



FIG. 9: Gaps Amin of spin-1 ladders and tubes with J = 1. 
Marks "SPA" and "QMC" mean "derived by the SPA" and 
"derived by the QMC simulation" respectively. The QMC 
data, which all achieve their thermodynamic-limit values, are 
calculated by Munehisa MatsumotO(fi^ The method of deter- 
mining the gap and the correlation length in the QMC simu- 
lation is, for example, explained in Refs. |5^ adnlGa. 
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4leg tu.bfe 
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FIG. 10: Spin correlation lengths in spin-1 ladders and tubes 
determined by the SPA or the QMC simulation. The QMC 
data is provided by Munehisa Matsumotoi— 



B. [0, 0] -field case 



less modes. [For example, on the c = 1 -I- 1 critical GS of 
two independent spin-i (or spin-1) AF chains with a uni- 
form field a rung coupling brings the hybridiza- 
tion of massless modes. As a result, a massless mode 
becomes gapful and only the other one remains being 
massless.] While, all magnon excitations have to possess 
a wave number kr in the (frustrated) tubes, due to the 
translational symmetry along the rung. Moreover, the 
tubes also have the reflection symmetry (Fig. . There- 
fore, we expect that when the lowest doubly degenerate 
magnons are condensed in frustrated tubes, these two 
symmetries strongly restrict possible interactions and hy- 
bridizations between the low-energy excitations. Conse- 
quently, it would be natural that two kinds of massless ex- 
citations are present after the lowest two magnon bands 
are condensed. It has been suggested in Ref. H^j that 
a strong uniform field brings a c ~ 1 + 1 phase in the 
3-leg AF-rung "spin-i" tube. (Note that this massless 
phase of the spin-i tube is located just under the satu- 
rated state. While, the massless phase in the integer-spin 
tube, expected here, appears just when an infinitesimal 
magnetization occurs.) Therefore, the emergence of a GS 
with two bosonic massless modes might be universal for 
odd-leg AF-rung (frustrated) spin tubes with any spin 
magnitudes. 

We investigate the above expectation of the c = 1 + 1 
phase in more detail, below. First, we demonstrate that 
our NLSM plus SPA method for ladders and tubes can 
reproduce the Zeeman splitting in Sec. HUB II Sub- 
sequently, in Sees. IIII B 21 and HUB 31 we discuss the 
magnon condensed state in frustrated tubes using a 
heuristic bosonization method. The discussion there 
takes into account the translational symmetry along the 
rung and the reflection one more carefully. The results 
further support the presence of the c = 1 + 1 state. 



In this subsection, we investigate the ladders and tubes 
with the uniform Zeeman term H2a|l . To this end, one first 
had better notice the following general aspects in ID spin 
systems, (i) The uniform field always splits triply degen- 
erate spin-1 magnon states into — 1, and —1 states, 
(ii) When a magnon band crosses the zero energy and 
the magnon condensation occurs in a ID U(l)-symmetric 
spin system, the GS is usually regarded as a c = 1 one- 
component TLLiiaSlSLSLSmii The combination of 
these statements and the band structures in Fig. 0] leads 
to the conjecture that as a sufficiently strong uniform 
field is applied, nonfrustrated systems enter in a stan- 
dard one-component TLL phase (c = 1), whereas the 
GSs of frustrated tubes become a two-component TLL 
(c = 1 + 1). This is another new even- odd property in 
AF-rung tubes. 

However, as multimagnon modes are condensed simul- 
taneously, interactions among the resulting massless ex- 
citations (TLLs) can be present in general. They could 
induce the hybridizations and gaps in a part of the mass- 



1. Zeeman splitting in ladders and tubes 



From Eqs. (|13|l and (|18|l . we can represent the Eu- 
clidean action of the [0,0]-field case, in the Fourier space, 
as follows: 



+J^^,A,K]+iNl3LX,p, (25) 



where the mark means Hermitian conjugate, J^aO<-) = 
(n? (k) , . . . , (k) ) , (k) = (^A^, (k) , ^A^^ (k) ) . The 
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2N X 2N matrix A±^ and N x N one are defined by in the action (|25|l . the effective action would be given as 



A^(k) 



A,(k) 



02 ai 



HuJn: 



(26a) 



■• a2 

a2 ai ) 



2gc^ 



where ai = (w,^j + e{k)'^)/{2gc), 0.2 = 



(26c) 
(26d) 

a2, ao = 

[02] for ladders [tubes], and 1 is an iV x iV unit ma- 
trix. We stress that A±_ is not Hermite but normal 
(A^Jl_l = A_lA^), so it can be diagonalized by a uni- 
tary matrix (hence, we do not need to consider the Ja- 
cobian generated from the diagonalization in the path- 
integral formalism). From Eqs. (|A2|) - (|A5|) . the eigenval- 
ues of A±_ [A^] are Af = (w^ + er{kf ± 2iHuJn)/{2gc) 
\Al = + €r{k)'^ + H'^)l{2gc)\. Using these, we can 
integrate out ni in Z and derive the SPE, 



27r ^ 



dk 



E 

=0,+, 



coth ( — e 



(^e^(fc)) = TV, (27) 



where e°(fc) = c^J B + ir"^ + m / and e±(fc) = eO(fc) T 

H. The quantities e°'+'~(fc) can be considered as the 
magnon dispersions. The comparison between Eqs. H27|) 
and H22(l elucidates the relation £,r{H ~ 0)^^ = 
S,r{H)~'^ + H"^ /c^ at T = 0. Thus, one can see that new 
dispersions e^(fc) restore the Zeeman splitting at T = 0. 
The magnons with e°(fc), (^) e~{k) take 5^ — 0, 
1 and —1, respectively. 

Strictly speaking, the above SPA cannot handle the 
strong-field case {H ^ J), where some magnons are con- 
densed, since the magnon condensation and the finite 
uniform magnetization are not taken into account in the 
Haldane mapping Q). However, as will be explained be- 
low, under certain approximations, the above considera- 
tion can be extended to the situation where the magnons 
are condensed. 



2. Magnon condensed state in frustrated tubes 

Referring the arguments in Refs. (gH and FFol we try 

to construct the effective theory for the magnon con- 
densed state in frustrated tubes. After integrating out 
the magnon fields except for the degenerate lowest ones 



9=±P 



1 

.2^0 



H ■ {niq X rflq 

gc 



V{mq)}, (28) 



(26b) where p = ^ {kp = ^n) and 



mr(x) = Urini{'x.) — \ — ^^sin(fcr' + 7r/4)n;, (29a) 

1=1 



V{mq 



2gc « 



^^{^^f^u\mq\\ (29b) 



In the potential Vijfiq), Ag — Amin (Amin is defined in 

the no- field case), c^^^ = —^2gc — ' ^^'^ introduced the 
biquadratic term M|77ig|'* to ensure the stability of the 
magnon condensed state [one can interpret that it origi- 
nates from the large- expansion for the 0{N) NLSM]. 
We perform the Ginzburg-Landau (GL) analysis (i.e., a 
mean- field theory) for the action H28|l . The minimum 
of the GL potential V{mq) is at rfiq = in the weak 
uniform-field regime H < Ag. On the other hand, as H 
exceeds Ag, the minimum is located in the field configu- 
ration, 



K) + K) 



a^2 



7J2 



4ugc 



A2 

« _ - 2 



TO^ = 0. (30) 



Namely, 77 = Ag is just the critical point between the 
condensed phase with a finite uniform magnetization and 
the noncondensed phase, within the GL analysis. For the 
condensed case iJ > Ag, let us introduce the new field 
parameterization , 



.q _„.g^.„.g - {rhq + m)e^^^'^. (31) 
Substituting Eq. (|31|) in the action H28|) . we obtain 



cff 



E 

q=±p 



1 
2^0 



+■ 



2gc 



gc 

■KW? 

2^ 



mqf + c^{d:,mqf 



ii'-^l -2 


gc 



(2mraq 4- m )drOq 

{dreqf + C^d^Oqf 



}, (32) 



up to the quadratic order of the fields. The action 
indicates that fields Wg and ffiq are massive, and it means 
that the low-energy limit of the magnon condensed state 
is described by the two phase fields 0g=±p. If one inte- 
grates out the massive fields neglecting the third or higher 
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derivatives of all fields in the action (|32(l . the resultant 
effective Euclidean Lagrangian density is 

K 
2^ 



E 

q=±p 



where K 
and /i2 = 



{3H 



— ihidrOq, 

Al)/{H^-Al). 



(33) 



hi = ^Hm'^l{gc), 
This is just the La- 



grangian density for a two-component TLLi^*^ K and v 
correspond to the TLL parameter and the Fermi veloc- 
ity, respectively. These two quantities are renormalized, 
from their values of the GL theory, by the interactions 
neglected here (we will discuss the effects of such inter- 
actions below). It has been known that in the vicinity of 
the lower or upper critical fields, the TLL parameter for a 
gapless AF spin chain preserving the z component of the 
total spin is generally approximated as unity correspond- 
ing to the free fermion.^^-^^-^-^ Under the assumption that 
this property also holds in tubes, K in Eq. is close 
to unity as — > Ag -I- 0. Here, imitating — Urini, 
we define the Fourier transformation of the spin operator 



as Tj- =^ Uj'iSi j 
and ifT^ serve 



Tr^n + (-l)^T,.,,t. Equations (EU, ®, 



-1= T7(^? X 



1 
1 



E—im 



■ - 2 

[1 771 



T+ 



q.st 



(34a) 
(34b) 



q is l/v47r in our notationia the equal-time 

q is defined as 
where 0= is the Heav- 



{mq -\- mje' 

where the latter equation suggests that the radius of 
Oq is I/a/tTj a-iid it also means that the radius of the 
dual field 

commutation relation between 9q and 
[(t>q{x),9q{y)] = iQsiy - x^, 
iside's step function. We expect that drOq [e^^^t) in 
Eq. H34|l is the most relevant bosonic term in the effective 
representation of Tq ^ (T^^J. However, it is difficult to 
determine the forms of second relevant or more irrelevant 
terms within only the above NLSM plus GL analysis. 

In order to examine a more proper bosonic represen- 
tation of Tj., let us once review the low-energy physics 
of a spin-1 AF chain with a uniform field. In addi- 
tion to the NLSM approach in Sec. ^ there is another 
low-energy effective theory for the spin-1 AF chains. 
The latter approximates the spin-1 chain without exter- 
nal fields as three copies of massive real fermions, each 
of which is equivalent to an off-critical transverse Ising 
chain In this scheme, the spin operator is writ- 
ten as 



s:/a 



S: + (-iys:„ (a = l,2,3orx,y,z), 

(a -I- 3 = a), (35) 



where £,2(r) the left (right) mover of the ath real 
fermion, Cct ifJ-a) stands for the order (disorder) field in 
the a-th Ising system, and Cst is a nonuniversal constant. 
The GS in the fermion picture corresponds to the disor- 
der phase in the Ising picture: (/ia) ^ 0. When the uni- 
form field exceeds the lower critical value, the low-energy 
physics is governed by a TLL with a scalar boson and 
its dual 9. In this case, the representation of the spin op- 
erator (|35|) should be changed by using these two fields. 
Equation (|35|l and the known results of the 2-leg spin-i 
AF ladder with a uniform fieldSSi^i provide a desirable 
representation. 



5,1 



Ci cos 



47r0 - 



27rMua; 



5. 



St 

a 



(72(73 cos 



ttM„ X 



a I 



-a 



4 sm 



2ttMuX 



, (36a) 
(36b) 

(36c) 



where — and Ci_4 are a constant. The 

third fermion system is still massive (it corresponds to 
(/is) 7^ 0) and the fermion ^ stands for the magnon 
with S*^ = 0. The formula ^ is vahd in M„ < l-ii 
Equation (|36c|) means that the radius of 9 is l/\/7r- 
While, Eq. Ij36b|l also suggests that the radius of cj) 
is the same value, IZ-v/tt- This apparently contradicts 
the framework of the standard TLL theory. Actu- 
ally, if we allow the presence of both e'^* and e*^"^, 
the commutation relation between these two is nonlo- 
cal: e'^'^'^^^e^^^*^^-' = Q-^l<^>{^),Siy)]^i^/1T9(y) ^i^4>{x) _ 

sgn(?/ - a;)e'^^(^)e'^'^(^) (a; ^ y), where sgn(y - a;) is 
the sign function. The nonlocal property is inconsistent 
with the fact that the original spin-1 operators are mutu- 
ally local (i.e., two spins on different sites commute with 
each other). However, observing Eq. (|36|l carefully, one 
can find that the nonlocality between 173 and /is and that 
between cos(Y^(/)-f-7rAf„x/a) and e'^^ cooperatively re- 
store the locality among the spin operators 5*1 and . 
Therefore, we believe that the formula H36|) is valid and 
the radius of 9 (0) may be defined as I/^/tt (I/VAtt). 
Here, further notice that the effective Hamiltonian in the 
spin-1 chain has to be constructed by the sum of terms 
being locally related with each other, because the original 
chain is a locally interacting system. This statement must 
hold in the effective theories of the ladders and tubes 

Now, we go b ack to the frustrated tubes. Following 
Refs.^ andrTOl one can see that the rfip (or m_p) part 
of the effective action (|28f) [or H32|l] is the same form 
as the effective one for the spin-1 chain with a uniform 
field. Moreover, the bosonic representation H34() is very 
similar to that of the spin-1 chain, Eq. H36|) . From these 
facts, it is expected that Eq. H36|l helps us improve the 
imperfect formula (|34|l . We propose the following new 
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bosonic representation of Tr'. 



.{d^(t>p + d^4>^p) + 2Mt 



-C(i|cos V^(l)p 



2Tr MfX 



rpZ 
9, St 

T+ 

9, St 



COS 



Ct2 COS 



27rAf,a; 



a 

irMtX 



(37a) 
(37b) 



C,3e'v^««|l 



-Cf4 sin 



27r Mtx 



}, (37c) 



where (j)q is the dual of 6q, and Cti-t4 are a constant. 
The first term in Eq. (|37a|l would be acceptable from 
the real-time operator identity ^dtOq ~ dx4>q- The pa- 
rameter Mt can be fixed by the magnetization per site 
{Si j) = 2Mt/N . (Determining the correct value of Mt 
is difficult within the GL theory.) To preserve the lo- 
cality among the spin operators , we should regard 
that Ct2 and Co contain massive fields such as and 
/i3 in Eq. H3t)|) . However, since we actually cannot de- 
termine what massive fields the constants Ct2,ts, contain 
within the present heuristic approach, Eq. (j37b|l might 
be somewhat doubtful. 

We proceed to the discussion employing the for- 
mula H37|l . So far, we have omitted the interactions gen- 
erated from the higher-order terms of fhq and the trace 
out of the massive fields. We hence study their effects 
towards the two-component TLL. The interactions can 
induce terms involving (j)q and 9q in the effective Hamil- 
tonian for the TLL. Let us concentrate on the situation 
near H — > Ag, and assume that K , the radius of 0g, and 
that of 6q are approximated as unity, 1 / \/^, and 1 / y^, 
respectively. In such a case, the relevant or marginal 



vertex operators are restricted to e**^ 
e±'^*±, and e±*^^"®± (n = 1 or 2), where we de- 
fined <i>± — (t)p ± (f)^p and Q± — 9p ± 9^p. [Note that 
the scaling dimension of a vertex operator g^'^^'f'qi^q) is 
A/4tt in the Lagrangian H33|) with if = 1 in our nota- 
tion.'*] It is sufficient to investigate whether these terms 
can be allowed or not in the low-energy effective theory, 
in order to know how critical state appears in the frus- 
trated tubesiS To this end, we utilize several symmetries 
in the spin tube systems. 

From Eq. H31|l . the U(l) spin rotation around the z axis 
corresponds to the transformation 9q ^ dq + constantiS 
Since the spin-tube Hamiltonian should be invariant un- 
der the U(l) rotation, the effective theory does not 
have any interaction terms with e^^V^"-^<i and e='=*^"®+ . 
Equation H37|) shows that the one-site translation along 
the chain is identified with (j)q ^ (j)q + \/n{Mt ±1) and 
9q ^ 9q± y/n.'^^ Thus, the appearance of e^^^^'^" and 
g±j%/47r$+ jg g^j^gQ prohibited as far as Mt is not equal to 
a special commensurate value. (For the nonfrustrated 



tubes or ladders, the restriction from the above two sym- 
metries is sufficient to confirm the c = 1 state.) 

To further restrict the possible terms of vertex oper- 
ators, we consider the symmetries with respect to the 
rung direction. All the tubes have the reflection sym- 
metry illustrated in Fig. 13 the corresponding transfor- 
mations are Sij Sn-lj, ni — > n^^i {Sn.j and 
are fixed), Tr — *■ T_r, and — > m_r. From Eq. I|37|l . 
these obviously require the effective theory to be invari- 
ant under the mapping {4>q,9q) {(j)-q,9-q). This pro- 
hibits the sine- type operators, sin(/30_) and sin(/3$_) 
[/? G M] , since they change their signs under the mapping. 
To discuss the remaining relevant terms cos(\/47r$_) 
and cos(-y/7rZi0_), we further examine the invariance un- 
der the translation along the rung: Si,.j — > Si^qj and 
ni fii+Q [I mod and Q € Z]. From Eq. H29a() and 
the definition of Tr , these operations cause 

Tr -> cos(Qfc^)T; - sin(Qfc^)f_^, (38a) 
rrir — > cos(Qfcr)mr — sin(Qfcr)m-r- (38b) 

Comparing Eqs. (|38() and H37() . we propose the following 
transformation for the vertex operators: 

g±iv^e, ^ cos(Q/cg)e±'^'''' -sin(QA:g)e±*^^-', (39a) 
g±iv^0, ^ cos(Q/fcg)e±^^'^'' -sin(QA:g)e=^*^'^-«.(39b) 

In fact, as far as one focuses on the most relevant term 
in Eqs. H37b|) and (|37c|) . Eq. is consistent with the 
transformation H38() . [We will discuss the second relevant 
term in Eq. I|37c|l later.] The transformation H39|l leads 
to 

e^'^^- ^ cos2(Qfcg)e±*^®- -sin2(Qfcg)eT'^®- 
+0±, (40) 

where 0± = 2 sin(2Qfcp)(e±*v^^+e=F»v^e)+ - 
giiv^e-gT'v^e-) jf negligible in the 

sense of the point-splitting technique, ^''^'^'^ Eq. 140(1 
provides 



cos(-v/7i'6-) — > cos(2Qfcp) cos(\/7r6-)- 



(41) 



Due to cos(2Qfcp) ^ 1, cos(-y7rB_) has to be absent in 
the effective Hamiltonian. Furthermore, using Eq. H40|l . 
one can obtain 

cos(2x/^e_) [cos''(Qfcp) +sin''(Qfcp)]cos(2V¥e_) 
-a, (42) 

where O = 2s\T?{Qkp)[{e^^^~ e'^^^- + 
g-iv/5Fe-giv^e_) _,_ (j^ ,,)] gy ^gjj^g ^Yie point split- 
ting, O may be replaced with a constant. Since kp can 
not satisfy cos'*(Qfcp) -I- sin''(QA:p) — 1 and sin(Qfcp) — 
simultaneously, the marginal terms cos(2y^O_) is also 
forbidden. From Eq. ((39b(l . the similar argument from 
Eq. (|40|l to H42|l . of course, can be adopted to the vertex 
operators with (j)q. 
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From these arguments, we can say that the symmetries 
of tubes make aU the relevant or marginal operators ab- 
sent in the effective theory. Namely, the above bosoniza- 
tion argument supports the presence of the c = 1 + 1 
state. 

Here, we had better think again the proposal (|39|l . The 
reader will immediately (or already) find that the final 
term in Eq. (|37c|) does not obey the desirable transforma- 
tion corresponding to Eq. H38|l . Therefore, it is expected 
that either the term in Eq. (|37c|) or the proposal Ij39b|) 
is invalid. In the latter case, one cannot forbid the exis- 
tence of cos(2-y7r$_). Then, instead of the discussion on 
the symmetries, let us count on the known result: in the 
single integer-spin-^ AF chain, the TLL parameter K 
monotonically increases together with the growth of the 
magnetization (5"!) within the region (S'J) < ^-.eSjTO^ 
Provided there exists the same nature in the integer-spin 
frustrated tubes, the scaling dimension of cos(2^/7r$_), 
2K , is larger than 2 (i.e., irrelevant) for the small- Af^ 
case. Therefore, we can predict again the presence of the 
c = 1 -f 1 state. 

Now, are there any vertex terms which survive from 
the restriction of symmetries of the reflection Tr — » T_j. 
and translation l|38|l ? We can find that the following four 
terms: 

cos(-\/47r0p) -t- cos(\/47r0_p), 
sin(\/47r'/'p) + sin(V47r0„p), 
cos(V47r6'j,) + cos(V47r6'_j,), 
sin(V4^6'p) + sin(V4^6'„p), (43) 

are invariant under these two operations. This is consis- 
tent with the presence of the final term in Eq. (|37a|l . 



3. Stronger uniform-field case in the frustrated tubes 
(quantum phase transition) 

We next discuss the frustrated tubes with a stronger 
uniform field, where the second lowest magnons are con- 
densed as well as the lowest ones. 

First, we consider the 3-leg tube. The second lowest 
magnon corresponds to the field mg. When its conden- 
sation takes place, the new phase field and its dual 0o 
would emerge from the field toq, like Eq. 131(1 . Because 
mo is invariant under the reflection n/ — > nAr_; and the 
translation along the rung l(38|l . the symmetries do not at 
all restrict the form of interaction terms with 6^ and 0o 
in the effective theory. On the other hand, other symme- 
tries of the U(l) rotation and the translation along the 
chain demand the invariance under ^ + constant 
and ((Ao, So) ^ (00 + V^(Mo ± 1), ± V^) Wo ^ Mt]- 
Therefore, all the vertex operators including 0o or Oq 
alone are prohibited. Are there any vertex operators with 
4>Q, 9q, (pg, and 9q which are invariant under all symmetry 
operations? Employing the term (|43|) . one can find the 



following terms permitted for all symmetries: 

cos[V47r(6'p - Oo)] + cos[V^{9^p - Oq)], 
sin[\/4^(0p - ^o)] + sin[V4^(0_p - ^o)]- (44) 

Notice that the same type of terms as Eq. (I44|l . where 
9q,o are replaced with (fiqfi, are not permitted because 
Mo 7^ Mt- Relying again on the known result that the 
TLL parameter K (Kq) for 9q {9o) are larger than (close 
to) unity, we can expect that the scaling dimension of 
terms (I44II . 1/K + I/Kq, is smaller than two, and they 
must be relevant. Introducing the new parameteriza- 
tion:^3 Go = {9p + 9o + 9,p)/V^, Gi = {9p - 9_p)/V2 
and 92 = {9p + 9^p-29q)/V&, one sees that the relevant 
terms ()44|l can be rewritten by the two fields 81.2 and 
do not contain the field 60- As a result, the field 60 
provides a one-component TLL, whereas the remaining 
two fields 8i_2 carry a gapful excitation. Thus, we can 
predict that the c = 1 + 1 state in the 3-leg tube are 
broken down to a c = 1 one once the condensate of the 
kr = magnon occurs. (Although the above new param- 
eterization makes the Gaussian part of three TLLs be a 
nondiagonal form, it would not influence the prediction 
of the c = 1 state.) 

The similar argument also holds in other frustrated 
{N > 5) tubes. In these cases, the second lowest bands 
are twice degenerate: the bands with the wave numbers 
fcp_i and /c_p+i. Thus, as the field is applied enough, 
two pairs of phase fields ('/'±(p-i), ^±(p-i)) appear cor- 
respondingly to the condensations of 'rn±(^p_iy Similarly 
to Eq. 1331), we can find the following interaction terms, 
which are permitted from all symmetry operations: 

cos[V47r(6'p - 6'p_i)] -I- cos[V^{9-p - 9p^i)] 
+ cos[\/47r(6lp - 9^p+i)] + cos[\/47r(6'_p - 9^p+i)], 
sm[VA^{9p - 6lp_i)] + sm[Vi^{9-p - 0p_i)] 
+ sm[VI^{9p - 0_p+i)] + sin[\/4^(0_p - 0-p+i)]. (45) 

These are expected to be relevant. As in the 3-leg 
case, if we introduce the new fields^ Qq = (6'p -|- 9-p + 
9p^i + 0-p+i)/\/4, 81 = {9p - 9^p + 9p^i - 0_p+i)/V4, 
82 = {9p - 9p-i)/V2 and 83 = {9-p - 0_p+i)/V2, the 
terms (|45|l are re-expressed by using only three fields 
©1,2, 3- Consequently, a c = 1 state with the scalar field 
©0 would appear, and other three fields have a massive 
spectrum. Thus, we can finally arrive at the general pre- 
diction that a c = 1 state emerges instead of the c = 1 + 1 
one as the second lowest bands crosses the zero-energy 
line in all the frustrated tubes. The quantum phase tran- 
sition between these two critical (c = 1 + 1 and c = 1) 
states would be observed as a cusp singularity in the 
uniform-field magnetization curve as in Fig. 1111 because 
the uniform susceptibility is generally proportional to the 
number of massless modes in ID spin systems. Moreover, 
the GS phase diagram for the frustrated tubes is drawn 
as in Fig. 1121 (Applying the same argument to nonfrus- 
trated systems, it is found that the c = 1 state lasts out 
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c=l phase 




Cusp singularity 



Uniform field H 



Gapped phase 

FIG. 11: Expected magnetization curves in frustrated integer- 
spin tubes. 




Ji 
FMrung AF rung 
(frustrated) 



FIG. 12: Expected GS phase diagram of odd-leg integer-spin 
tubes around {J±,H) — (0,0). 



even when the condensate of the second lowest magnon 
occurs.) Besides the above scenario of the magnetization 
cusp, other cusp singularities have aheady been found in 
theoretical studies of ID quantum systems i?-^i?^i??i?"^i?^'^-?^ 
It is known that around such cusp points, the left or right 
derivatives of the magnetization (i.e., susceptibilities) al- 
ways diverge. However, such a singular phenomena is 
expected not to occur in our cusp mechanism. Thus, we 
may insist that the cusp in the frustrated tubes is a new 
type. 

Finally, we remark the validity and the stability of the 
bosonization arguments in Sees. 1111 B 21 and |lll B 31 The 
proposal (|39|l is a strange form: other symmetry opera- 
tions correspond to a transformation for the phase fields 
own, but Eq. (|39|l is the transformation for the vertex 
operators. In fact, as stated already, Eq. H39I) is incom- 
patible with the final term in Eq. (|37c|) . There might be 
a more natural transformation than Eq. (|39|l . As far as 
we know, it has never been elucidated whether the GL 
theory is efficient or not even when multimagnon bands 
are condensed. Furthermore, the coupling constants of 
the phase-fields interactions [for example, Eqs. 1)44(1 . (|45|l . 
and other irrelevant interactions with phase fields] might 
be as large as the order of J. In such a case, the pertur- 
bative treatment of them and the biquadratic term u|m|* 
would be dangerous. Thus, there is still a little possibil- 
ity that the interactions break down the c = 1 -I- 1 state. 



On the other hand, we essentially use only four symme- 
tries of the spin tubes in order to lead to the c = 1 -I- 1 
state. Thus, if our strategy based on the bosonization 
and the GL theory is admitted, the c = 1 -I- 1 phase must 
be stabilized against several small perturbations preserv- 
ing those symmetries (e.g., XXZ-type anisotropy, single- 
ion anisotropy -D ^■(5'f j)^, next-nearest-neighbor cou- 
pling, etc.). 



C. [0, 7r]-field case 

This subsection discusses ladders and tubes with the 
[0,7r]-field term l|2b|l . Odd-leg tubes, which include frus- 
trated tubes, are not admitted in this case. 

Imitating the argument in Sec. 1111 B II 
also express the action with the quadratic 



bosons 



diag(-l,l. 



which 
■■•,(-1)^) 



we can 
form of 
F — 

In the action, the 2N x 2N 



is replaced with 



gc 

matrix held between and A/j^ are not normal due to 
Fi^^ . Its diagonalization hence generates a nontrivial Ja- 
cobian differently from the [0, 0]-field case. To avoid this 
difficulty, we turn to the real-time formalism, even though 
it restricts our consideration to the zero-temperature 
case. The partition function Z = J Vni'DXe^^E is as- 
sociated with the real-time vacuum-vacuum amplitude 



Zv = / VniVXe 



via S'p''^'[t] — iS-^ 



[0,7r] 



it]. In 



the Fourier space (note that the frequency w is a real 
number, and not the Matsubara one w„), the real-time 

action S^^'^^ can be represented as 



4°-^'[{n,},A,p] 



E 

k 



NTLX^p, (46) 



where k = {uj,k), Af^ = «(k), 



'AT 



(k)), Aff' = 



(^Afl^, ^M^), T is the time distance from the initial vac- 
uum to the final one, and we performed a unitary trans- 
formation rv^ — > (— l)'n]' for convenience. Two matrices 



and are defined as 



A«(k) 



— E,,, Br, 



^k-:^l, 

25c 

Huj- 
-I 1, 



(47a) 
(47b) 
(47c) 



where A^ [B^ is the same form as Ak, in which fii — s- 
(w^-c^/c^-f 2gcAsp)/(2gc) and ao,2 -ao,2 [ao,2 -> ao,2]- 
Since A^ and A^ are Hermite and can be diagonalized 
by a unitary matrix, we can apply the SPA prescription 
as in the preceding subsections. Eigenvalues of A J and 
Af are, respectively, Af = (w^ - e{kf)/{2gc) ± /^(k) 



and = (w^ — er(fc)^ — H^)/{2gc), where we defined 
e-(fc) = c^k^+l-\ ^--2 = -2.gAsp/c, 
= c^k^+ir\ = r' + §^ cos fc„ 

I 



where 



Here, we used the so-called ie-prescription^ in the 
derivation of the SPE (gHl- The SPE indicates that Asp 
is real (not imaginary) and negative in contrast to the 
imaginary-time schemes. It is verified that as iJ — > 
[Jj_ ^ 0], Eq. gni) is reduced to the SPE ^ [d] of 
the zero-temperature case, where (^, Asp) corresponds to 
(?,«Asp). 

One can regard e^'^(fc) as the magnon band disper- 
sions. At iif = 0, ^{k), ^ri^): ^^'^ coincide with 
tr{k) in Eq. (|23l, e|r|(fc) and eN+i~\r\{k) [e[N]+i_|r| (fc)], 
respectively, for ladders [tubes: > 3]. Although it is 
hard to solve the transcendental equation H49|l for Asp, 
we can extract some features of the band splitting in- 
duced by the [0, tt] field, from Eqs. (03 and (|S3). The 
inequalities e^(fc) > e~(fc) and hr{k) > 2H^ > show 
that the final term in the left-hand side of Eq. (|49|l is 
negative or zero. Therefore, decreases and the bands 
e~(fc) fall down as H is applied. The form of the dis- 
persion (|50a|l indicates that as H is applied, the triply 
degenerate bands er{k) with cosfcr > (< 0) are split 
into doubly degenerate upper (lower) bands and a non- 
degenerate lower (upper) one. The former two are the 
transverse modes, and the latter is the longitudinal one. 
The bands with /c^ = ±-|, which are present only in odd- 
leg ladders and (4 x Q)-lcg tubes [Q £ Z], are divided 
into three bands. From these considerations, we can il- 
lustrate the band splitting as in Fig. ^| [Remember 
that tubes have a sixfold degeneracy in the no-field case 
(Fig. 0J.] The figure tells us three remarkable aspects, 
(i) Any band 6^!'+'" split by the [0, tt] field tends not to 
approach the other neighboring bands (to avoid the band 
crossings). This contrasts with the Zeeman splitting in 
the uniform ([0, 0])-field case, where the crossings among 
magnon bands with different indices kr occur, (ii) The 
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Employing these eigenvalues, one can trace out fii in Z^, 
and then obtain the effective action S^'^^ [Asp]- The SPE 

54°'"'^p]/aAsp = o can be calculated as 



N, (49) 



(50a) 
(50b) 




FIG. 13: Predicted band splitting in the [0, 7r]-field case. 



lowest bands are doubly degenerate in all systems. It 
might imply that a sufficiently strong [0, tt] field always 
engenders a c = 1 -I- 1 phase, (iii) The FM-rung coupling 
competes with the [0,7r] field. However, the figure sug- 
gests that any qualitative differences between the com- 
petitive and noncompetitive cases do not emerge at least 
in the weak rung-coupling regime. 

We believe that the band structure in Fig.^|is qualita- 
tively valid. However, its details would strongly depend 
upon the SPA strategy. Particularly, we mind that the 
symmetry corresponding to the remaining double degen- 
eracy of transverse modes cannot be found. Therefore, 
the degeneracy and the prediction (ii) would be ruined 
in more quantitative approaches. 



D. [7r,0]- and [tt, 7r]-field cases 

This subsection addresses the [tt, 0]- and [tt, tt] -field 
cases, in which the external field is alternated along the 
chain. Utilizing Eqs. ljBl|) and H18|) . we can describe their 



277 ^ 



dk 



2H^ ( 1 1 



,X^^_~^r{^) hr{k)Wik) Irik) 



{k) = ^Jlr{kY + m, ef{k) = ^e{kY + 2m±hr{k), 



hr{k) = J 4H^ + 4H^e{ky 



( 2J^ 
V Ja2 



cos kr 
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low-energy action as follows: 



SE[{fii},X] 



dx 



I ■ ni 



-iiVA(x) 



(51) 



where Ki H/a {Ki = {-lY+^H/a) for the [tt, 0]- 
field ([7r,7r]-field) case [H = {0,0, H)]. Notice that 
odd-leg tubes are not permitted for the [tt, 7r]-field case. 
From Eq. 1)51(1 . it is found that a unitary transforma- 
tion n/=cvcn — —ni=cvcn exchanges the action of the AF 
(FM)-rung system with the [tt, tt] field into that of the 



FM (AF)-rung one with the [tt, 0] field. Therefore, it is 
enough to investigate only either [tt, 0]-field case or the 
[tt, 7r]-field case. We analyze the former case below. 

After diagonalizing the quadratic part of ni in Eq. 151|l . 
the action are arranged as 



S Jr ' 



.| -Ki7VA(x) 



(52) 



where 



Ar(x) 



1 r 



-2i.gA(x) 2Jj_ 
c Ja^ 



cos kr 



2{N+V, 



-1 - 
H 



(ladders) 
(tubes : iV > 3), 



(53a) 
(53b) 



and the field mr(x) = Urini [see Eq. (|29a|l ]. The action 
(|52|l can be considered as a generalization of that of the 
single chain with a staggered field (|B1|I . Thus, we can 
apply the Green's function method in Appendix IbI to the 
present [tt, 0]-field case. (We would like the reader to 
refer Appendix IbI or Ref. (33 before proceeding below.) 
Following Appendix IbI let us define Green's functions: 



A,(x)GO(x-x') 
G'^(x-x') 
G'^(x-x') 



52(x-x'), (54a) 
(r,TO^(x)m^(x'))c, (54c) 



where the subscript c means "connected," 7^ denotes 
imaginary-time ordered product [see Eq. ((B6p ]. and 
5^(x — x') = 5{x — x')5{t — t'). In anticipation of remov- 
ing the space-time dependence of S^r via the SPA process, 
we already assumed that the above Green's functions de- 
pend only on the distance between x and x'. The magnon 
dispersions of transverse and longitudinal modes can be 



As in Eq. (|B10|) . the Fourier transformation of is es- 
timated as 



G°(k) 



(55) 



COS fcr and Asp is the saddle- 



where = 

point value of A(x). Using and referring the way 
deriving Eqs. (|B4(I and l|Blip . we obtain the following 
SPE determining Asp and <^r, 



TV- 



(56) 



where the final term in the right-hand side represents 



the [ 7r,0' 
Eq. (IHTTll 



-field effect. As J±^ 0, Eq. 1(56(1 returns to 
Applying Eqs. I(B5() and ((B12() . we further 



determined from G^(x— x') and G^(x— x'), respectively. evaluate the staggered magnetization as 



Ml 



S{ni) = Uriirfir 



N+i Sr=i '^r,odd sin ^ ^'^^ ^ tan 2( w+i) ) 



2 - 
H 

J 



H 
J 



(ladders) 
(tubes : iV > 3) 



(57) 



where M/ is parallel to the field H like the single-chain case [see Eq. ((B12ll ]. namely M; = (0,0, M/). The stag- 
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0.8h 4,5,6,7-leg 




0.1 0.2 

H [7t,0] field 



FIG. 14: Staggered magnetizations of A^'-leg spin-1 FM- or 
AF-rung tubes with the [tt, 0] field, A'/f The dotted curve 
A/chain stauds for the staggered magnetization of the sin- 
gle AF chain with the staggered field (see Fig. 0. The 2- 
leg tube means the 2-leg ladder. A relation Mf-'°s(i^) < 
^^(p+i)-ieg(.j:^) [MP-'°s(//) > MiP+''"''=^(//)] is realized for the 
FM [AF]-rung case. Magnetizations Mf'^''^''^'^"^ almost over- 
lap. 



gered magnetizations are independent of the chain index 
I in the tubes: Mi = Ms- On the other hand, the l- 
dependence of the magnetizations clearly exists in the 
ladders. We emphasize that these inhomogeneous dis- 
tribution of the staggered magnetization cannot be pre- 
dicted by the standard NLSM scheme, which assumes 
a short-range AF or FM order to arise for the rung di- 
rection. From sin(j^) = sin C'^^+j^'^^ ), we find that 
Ml = Mn+1-i is realized in the ladders. The results 
Ml = Ms in tubes and Mi = M^+i-i in ladders indicate 
that the present NLSM plus SPA scheme preserves the 
translational symmetry along the rung in tubes, and the 
reflection one about the plane containing the central axis 
of the ladder, in the [tt, 0]-field case [refer to the discus- 
sion about the validity of Eqs. lfTB|) and lfT7|) in Sec. 1111 Aj . 

We will investigate the magnetizations Mi and the 
magnon modes in detail, below. 



1. Staggered Magnetizations 



If Asp is fixed by the SPE, the magnetizations Mi are 
also done in Eq. H57|) . Figures Fl 41 1 61 displav them in zero 
temperature. The FM rung coupling and the [tt, 0] field 
are cooperative each other, while the AF rung coupling 



competes with the field. A consequence is, as shown in 
Fig. El that for FM (AF)-rung tubes, Ms are always 
larger (smaller) than that of the single AF chain with 
a staggered field. In addition, the growth of N grad- 
ually enhances (reduces) the magnetization Ms in FM 
(AF)-rung tubes. The magnetization profile such as the 
panel (a) in Fig. 1141 is also expected in nonfrustrated 
(FM-rung) ladders. Actually, as expected, Fig. ^1 indi- 
cates that for the small-field regime H < 0.05, Mi tend 
to increase together the growth of N. It further explains 
that the more the Ith chain approaches the center of the 
ladder, the larger its magnetization Mi becomes. This is 
understood from the consideration that the chains near 
the center are more subject to the FM-rung correlation 
effects than ones near the edge. On the other hand, one 
can extract following two unexpected features in Fig. 1151 
(i) The maximum magnetization in Mi of the A^-leg lad- 
der, M^^^, is slightly larger than that of {N + l)-leg 
ladder for the regime H > 0.05 and > 3. (ii) The edge 
magnetization Mi^n is smaller than that of the single 
chain in the same regime. Since the [tt, 0] field and FM 
rung coupling, which is absent in the single chains, must 
cooperatively enhance the growth of Mi, these two re- 
sults are expected to be incorrect. Because in the spin- 1 
AF chain, the staggered magnetization obtained by the 
SPA is almost consistent with the DMRG data within 
< -J < 0.5,"^^ these unexpected results would mainly 
originate from the averaging of constraints, Eq. (|17() . The 
approximation (|17|l would prevent AIi from increasing in 
the regime H > 0.05. The true magnetization curves 
are expected to have a less Z-dependence so that the edge 
magnetization Mi^m is always larger than that of the sin- 
gle chain, Mehain in Fig. CHI Moreover, M^+i > M^i|^^ 
must hold for all region < iJ < oo in the FM-rung lad- 
ders. 

Figure lTBI provides the staggered magnetizations in the 
frustrated AF-rung ladders with the [tt, 0] field. The left 
panel [a] insists that the staggered magnetizations tend 
to point to the same direction as the [tt, 0] field, when 
the AF rung couphng is small enough; ^ < 0.03. The 
edge magnetization Mi jv increases most rapidly in such 
a weak rung-coupling region since the edge chain receives 
the competitive AF rung coupling from only one side, un- 
like other chains. The rapid growth of Mi,m and the AF 
rung coupling would make the growths of M/=even slower. 
While, the panel [b] contains the following interesting 
phenomena: when the [tt, 0] field and the rung coupling 
are sufficiently small and large, respectively (-j < 0.1 
and ^ > 0.03) in the odd-leg tubes, the field induces 
the staggered magnetization pointing to the opposite di- 
rection to it in the even-Z chains. This result is unique for 
the ladders, and does not appear in AF-rung tubes (see 
Fig. I14|l . Such a magnetization configuration staggered 
along the rung does not also occur in even- leg ladders, 
because the configuration cannot be compatible with the 
edge magnetization turning to the [tt, 0] field. One may 
call the result in the panel [b] as an even-odd property in 
the ladders with the [tt, 0] field. Although the panel [b] 
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spin-1 FM-rung ladders (J,JL)=(lr0.05) 
0.8 




0.2 0.3 0.1 
H [7t,0] field 



0.3 



FIG. 15: Staggered magnetizations of A'^-leg spin-1 FM-rung 
ladders with the [tt,0] field and J±/,J = -0.05. For I < 
[f ], Mi+i{H) > Mi{H) is realized. Similarly to Fig. d the 
symbol A/chain means m^. 



further implies the simultaneous crossings of (MijMa) 
and (M2, M4, zero- magnetization line), they might be a 
coincidence depending on the approximation H17|) . From 
the discussion about Fig. ^1 Fig. ^| might be also less 
accurate for the regime H > 0.05. 



2. Magnon Dispersions 

We next study the magnon dispersions of the [tt, 0]- 
field case. Following the functional derivative technique 



such as Eq. (|B6|I . we can represent Gr 



T{L) 



as 



G^(^)(x-x') = 5'Gj!(x-x') + 252y dy 
GO(x-y)GO(y-z)J,^(^)(z)(.-iJl^) 



(58) 



In the calculation of 5X/6J", we suppose that each J" 
is an independent external field. The relation J^'^ — 
leads to G^(x) — 5^G°(x). The Fourier transformation 
of is 



G^(k) 



S'gc 



(59) 



Therefore, we can immediately conclude that the disper- 
sion of the r-th transverse mode e^(fc) is given by er{k). 
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FIG. 16: Staggered magnetizations Mi of A'^-leg spin-1 AF- 
rung ladders with the [tt, 0] field. 



Similarly to the single-chain case, this mode has the dou- 
ble degeneracy corresponding to the x and y components. 
Furthermore, like the no-field case, the tubes exhibit the 
four-fold degeneracy Cr = e_r except for r = and 
modes. The SPE Q tells us that = -H!^ is 
enhanced by the field H (or J^). The transverse bands 
ej{k) and gaps = (0) hence rise monotonically 
with H increasing. 

Although the estimation of the longitudinal mode 
e^(fc) is rather complicated due to the presence of J,^, 
it is possible through the application of the method of 
deriving eL(fc) in the single chain with a staggered field. 
X as'"'"' 

The trivial relation -gjri — — ) = is available as the 
integral equation determining SX/6J^ in Eq. (|58|l . Imi- 
tating Eqs. (|B7p - ljB9|l . we can transform it as follows: 



J (iy/(x-y)(^i 



'^A(y) 



= -2M^G^(x-x'), (60) 



where 



Mr = S{ml) = 52G°(0)ff/a, (61a) 
/(x-y) = ^[6F,(x-y) 

r 

+ {Mr/SfGl{^-y)\, (61b) 
F,(x-y) = G^(x-y)G5!(y-x). (61c) 

Equations H58|) and (|60|) lead to the following expression 
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for the Fourier transformation of G^: 

E,|3r,(k) + 2A<JG,(k)| 

where Mr — Mr/S = {111^). The longitudinal mode 
ef:{k) can be fixed by the pole structures of the real- 
time retarded Green's function G^(k) = G^(k)|i^^^_iz, 
where k = (z, fc), z = uj + irj, and 77 — > +0. Here, as in 
Eq. I|B15|I . let us introduce new symbols, 



Gr{z) 


= G^'(k)/(5V), 


(63a) 


Friz) 


= 3fp(k)/(25c), 


(63b) 


rtot(z) 


r 


(63c) 




= m.c{3m)T{z), 


(63d) 



where we omit the subscript fc, and J-{z) is an arbitrary 
function of k. In terms of these symbols, we obtain the 
simplified expression of the real-time Green's function. 



Griz) 



Mi 



(fep 



(e.(fc)2-z2)[r,„,(,) + ^^_^ 



(64) 



Analyzing Eq. (|64|l . one can find the longitudinal magnon 
dispersions. 

[Tubes] The calculation of Gr[z) in the tubes is easier 
than that in the ladders, because the tubes take cx (5^,0 
and M-r oc (5r-,o- These properties bring 



Gr=^Q{z) = 



1 



AW 



cx G^(k). 



(65) 



The longitudinal dispersion e^(fc), thus, is equivalent to 
the transverse one e^(fc) for r 7^ 0. Namely, in the tubes 
with the [tt, 0] field, the rth magnon mode is triply de- 
generate like the zero- field case (Fig. QJ, except for the 
0th mode. Of course, there exists the additional degener- 
acy tr — c-r- On the other hand, the 0th mode Green's 
function is written as 



Go(z) 



rtot(z) 



(eo(fc)2-z2)rt„t(z)+X2 



(66) 



The form of Gq(z) is quite similar to that of G(z) in 
Eq. ||BT7|| . Moreover, at T = 0, rj'2(z) have the same 
form as F^'^ fixed by Eq. ljB16p . Therefore, following the 
calculation from Eq. (|ijl8() to Eq. (|ij22() . we can achieve, 
at T = 0, 

e^ikf = ,l[kf+MllV\,,{e]^{k),k), (67) 

where we restore the subscript k in . Note that in the 
derivation of Eq. l(H7|l. we assume (^(k) < 2e,nin(^/2), 
where emin(^) is defined as the minimum of all the trans- 
verse dispersions er{k) (see Fig. ^J. Because the in- 
equahty Tl^^{eQ (k) , k) > is realized under eg (A:) < 



FM-rung tubes 




0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 



AF-rung tubes 




0.1 0.2 0.3 



0.2 0.3 0.1 0.2 0.3 
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FIG. 17: Transverse and longitudinal gaps of A^'-leg spin-1 
tubes {N > 3) with the [■7r,0] field and J = 1. As a compari- 
son, we also draw the gaps of the spin-1 chain with staggered 
field (see Fig. |5J in the left upper panel. 



(|67|l explains that the 0th longitudi- 
is always larger than the transverse 
one eo{k). Employing the exphcit form of F^(eQ(0),0) 



2e„,i„(fc/2), Eq 
nal band 



see Eq. (|B2ip ]. we can calculate the longitudinal gap 
£0(0) as follows: 



^0 



^0 



where 



2Ml 
3g 



Ar^i^(A^/A^) 



K{x) 



2t:x^\ - a;V4 



: arctan 



2^/1 -a;V4 



(68) 



,(69) 



and we, of course, assumed that the minimum point of 
the band e^(fc) is located at fc = 0. The SPE and 
Eqs. H65|l - H68|l enable us to know all the magnon band 
structures in the tubes (A^ > 3) with the [tt, 0] field at 
T = 0. Figure ini shows the gaps A^'^ at T = 0. In this 
figure, we see that A^ — = A^ holds for r 7^ 0. 
Since the strong rung coupling destroys the condition 
eQ(fc) < 2emin(fc/2), our scope in Fig. El is restricted 
to the extremely weak rung-coupling regime. We, how- 
ever, believe that the gap behavior in Fig. 1171 is robust 
even with a moderately strong rung coupling. The lowest 
(highest) band eg is split by the field in the FM (AF)-rung 
tubes. Thus, there are the magnon-band crossings only 
in the FM-rung tubes. The manner of the lowest-band 
splitting goes for that of the single chain with the stag- 
gered field (see the upper panels in Fig. El E^^id Fig.j^J. It 
implies that an A^-leg spin-S* FM-rung tube has the same 
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low-energy properties as the spin-(iV x S) single chain 
even for the weak rung-coupling regime. The growths of 
gaps in the AF-rung case are slightly slower than those in 
the FM-rung case. It must reflect the frustration between 
the rung coupling and the field. As already predicted, one 
can verify from Fig. El that all the gaps monotonically 
grow up with H increasing. We, thus, conclude that in 
tubes the [tt, 0] field induces no critical phenomena at 
least in the weak rung-coupling regime, irrespective of 
the presence of the frustration. Although the SPA pre- 
dicts that only the degeneracy of the bands Eq is lifted 
by the field, actually the other bands are also expected 
to more or less split (because any mechanisms preserving 
the triple degeneracy of magnon bands are not found in 
the [tt, 0]-field case). 

[2-leg ladder] Let us next investigate the longitudinal 
magnons for the ladders. For the 2-leg case, the estima- 
tion of Gr is as simple as that in the tubes because of 
J| = 7VI2 — 0, which leads to 6*2(2) oc 6*^(6;). Therefore, 
the dispersion of the longitudinal mode £2 (fc) is identical 
with e^ik). On the other hand, Gi is written as 



Gi(z) 



Ftotl^) 



(70) 



Because the form of Gi is same as Eq. (|66|) . the process 
from Eq. (|66|l to Eq. (|69l) can be directly applied as a 
way determining the longitudinal dispersion ef'(fc). As a 
result, at T = 0, we obtain two equations. 



(71a) 



2Ml 



2 

E 



A 



■K{Ai/A^) ,(71b) 



under the condition ef (fc) < 2emin(fc/2). Equation H71a() 
indicates [2emin(fc/2) >] ef'(fc) > (fc). 

[3-leg ladder] Like G2 in the 2-leg ladder, M2 = 
leads to £2 (fc) = e|'(A;) in the 3-leg ladder. On the other 
hand, Gi_3 are more complicated than the Green's func- 
tions in the 2-leg case. After a simple calculation, they 
are represented as 



A{uj) + iB{lo) ' 



(72) 



where 



B 



E,E^Vl,,{z) - 2t^uj{E, + E^)Tl,(z) 



Cl(3) 
^1(3) 

Er 



-2riuj{Ml "1" ^^'3 
E,^,^Tl,iz) + 2r,u;Tl,{z)+Ml(^, 
E3^r)rldz)~2r,u;Tl,iz), 



(73b) 
(73c) 
(73d) 
(73e) 



Under the condition 



< 2emin(fc/2) and T = 0, in 
= 0, G1.3 are fairly simplified. 
In that case, the explicit forms of their imaginary part 
are 



which Tl, > and F^^^ 



GUz) = 



2riuj 



{Esrl.+Msr+MlMl+Si 



{E,+E3)rl,+Ml+Mi 



(74) 



where £1 = Ar|'^LL}'^{^l^^)'^ and £3 = 4ry^w^. At the limit 77 — > -1-0 [see Eq. (|B19|) ]. both Gl and G§ take the same pole 
structure as follows: 



lim 3 cx 5(g(w)), g^ui) 



EiEsTlM + MjEs + MjEi 
iE,+E3)rUu;)+Ml+Ml' 



(75) 



From the solution of g{ijj) = 0, we obtain the following two longitudinal bands e^(fc): 



4(fc)' = I Wikf + e^ikf] + :^{Ml + Ml) 

^ '- J Zi tot 



± 



1 



tot 



ie^{kr + esikr)Tl,+Ml+M 



ei(fc)2e3(fc)2ri„t + A^2g3(fc)2 + A^lei(fc)^ 



(76) 



where Tl^^^ means rl^^{e^{k), k). From this result, it is clear that at least within the SPA scheme, the [7r,0] field 
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engenders the hybridization between two magnon bands 
ei and £3 in the 3-lcg ladders. Provided that the mini- 
mum of (fc) is in fc = 0, the gap = (0) can be de- 
termined by the replacement (e^, ei, ea) (A^, Ai, A3) 
in Eq. H76() . In the zero- field limit H ^ 0, where 
Mi,3 0, Eq. ^ reduces to 



FM-rung ladders 



ei(fc)' + e3(fc)'±ki(fc)'-e3(fc)'| -(77) 



This result and the inequality ei(fc) < e3(fc) [ei(fc) > 
e3(fc)] in the FM [AF]-rung ladder reveal that and 
are, respectively, split from £3 and ei (ei and €3) in the 
FM (AF)-rung ladder. Therefore, and should be 
rewritten as £3 and ef (ef and £3 ) for the FM (AF)-rung 
ladder. The gaps A^ also can be redefined. 

[4-leg and higher-leg ladders] The logic calculating 
the longitudinal dispersions in the 2- and 3-leg ladders is 
successful even for the 4-leg ladders. We mention only 
the results of them. The identities A42,4 — leads to 
^2 4:{^) — ^24{^'^)- While, 6*1^3(2) take the same form as 
Eq. ifT^ under the condition |u;| < 2£niin(fc/2), except 
that FtQt — J2l=i and kr — ^ are replaced with 
J2t=i^r and kr — respectively. Therefore, ei^^ik) 
can be fixed like Eq. i|76|l . As easily expected, the evalu- 
ations of the magnon dispersions in the higher-leg ladders 
demand the more complicated analyses. Here we do not 
perform them. In principle, one can study all the longi- 
tudinal bands using Eq. 

We show the gaps of 2, 3, and 4-leg ladders in Fig. 1181 
The gap behavior is much similar to that of the tubes in 
Fig. El the highest (lowest) band is largely split into the 
doubly degenerate transverse bands and the single longi- 
tudinal band in the FM (AF)-rung ladders. The splitting 
of £3(fc) is considerably small. Namely, A3 and A^ al- 
most overlap. Observing carefully the numerical data of 
3-leg and 4-leg FM-rung ladders, we see that A3 (A^^) 
is a little larger than Af' (Af) for the case Af < A3 
(Af > A^"). Even though the band crossings in FM-rung 
tubes (Fig. I17|l are allowed from the translational sym- 
metry along the rung, the ladders do not possess such a 
symmetry. Therefore, the level crossing in Fig. 1181 might, 
in fact, be an avoided crossing. Moreover, more quanti- 
tative analyses would lift the remaining triple degeneracy 
of the bands £1^2 in Fig. ^1 

Summarizing all the discussions about the magnon dis- 
persions, we can conclude that the [7r,0] field engenders 
the monotonic raise of all magnon bands, and cannot in- 
duce any critical phenomena at least for the weak rung- 
coupling regime. While, we have already predicted that 
the other staggered field, the [0, tt] field, induces the gap 
reduction fFig. I13|l. Therefore, our results in the [0,7r], 
[7r,0], and [7r,7r] -field cases suggest that the spatial di- 
rection with the staggered component of the fields es- 
sentially affects the low-energy excitations of the spin 
ladders and tubes. 

Finally, we notice again that through the transforma- 
tion n;=cvon — ^ ~ni=cvcn, 

all the resuhs in Figs. I14I18I 
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FIG. 18: Transverse and longitudinal gaps of TV-leg spin-1 
ladders with the [tt, 0] field and J = 1. The gaps A^^ and A3 
almost overlap. 



can be interpreted as those of the [tt, 7r]-field case. In ad- 
dition, note again that odd-leg tubes are absent in the 
[tt, 7r]-field case. The [tt, tt] field competes with the FM 
rung coupling. This frustration, of course, can induce 
the even-odd property as in Fig. El The tubes with a 
[tt, tt] field do not have the one-site translational symme- 
try along the rung, and do only the two-site one. There- 
fore, the band degeneracy caused from the translational 
symmetry should partially vanish in the [tt, 7r]-field case. 
However, due to the symmetry restoration via the map- 
ping n;=ovcn —ni^even, which would be valid only in 
the low-energy limit, the tubes with a [tt, tt] field take the 
same bands as those of the tubes with a [0, tt] field within 
our strategy. 



IV. SUMMARY AND DISCUSSIONS 

We provided a systematic analysis for the low-energy 
properties of A^-leg integer-spin ladders and tubes 
with several kinds of external fields (|2Jl within the NLSM 
and SPA framework. Our results would be reliable for 
the weak rung-coupling, weak external-field, and small- 
N cases. Furthermore, we expect that several results are 
robust even for the strong rung-coupling, strong external- 
field, and large- Af cases. Although we concentrated on 
only the zero-temperature case, the SPA strategy used 
here, of course, can be applied to the low-temperature 
case. 

Our results are summarized as follows, (i) In the 
no-field case, we derived the magnon band structure in 
Fig.^ and predicted a new even-odd nature: for AF-rung 
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tubes, only the odd-leg tubes possess the sixfold degener- 
ate magnon band as the lowest one. The sixfold degener- 
acy is not a merely approximate result, and is protected 
by the translational symmetry along the rung. Several 
SPA results were compared with the QMC data in Figs. El 
and 1101 (ii) In the [0,0]-field case, we predicted another 
even-odd nature: when the field is sufficiently strong and 
a finite uniform magnetization emerges, the GS of odd- 
leg AF-rung tubes becomes a c = 1 + 1 massless state 
(two-component TLL), while a standard TLL state with 
c = 1 appears in other systems. Generally, Zamolod- 
chikov c theorem^ prefers the emergence of a c = 1 state 
to that of higher-c ones in ID U(l)-symmetric systems. 
However, we predicted, using the GL and bosonization 
analyses, that the translational symmetry along the rung 
and the reflection one (Fig. ^ in the frustrated tubes 
make the c = 1 + 1 state stabilized. Inversely, once these 
symmetries are broken down (for example, due to an in- 
homogeneous rung coupling), the c = 1 -I- 1 state would 
disappear and a. c — 1 state emerges instead. Regard- 
ing the case where the uniform field is further strong, we 
also predicted that the above c 1 + 1 state is taken 
over by a c = 1 one when the second lowest magnons are 
condensed. At the transition from the c = 1 -I- 1 state 
to the c = 1 one, one could see a new cusp structure, 
which does not accompany the divergence of the suscep- 
tibilities, in the magnetization curve. Furthermore, the 
validity of our GL theory was briefly discussed, (iii) In 
the [0, 7r]-field case, the SPA analysis suggested that the 
lowest doubly degenerate bands go down with the field 
increasing, in all systems. From this, one may think that 
a c = 1 -I- 1 state is also possible in the [0, 7r]-field case. 
However, it is doubtful since we were not able to find 
any symmetries leading to the degeneracy of the lowest 
two bands. We thus anticipate that the above double 
degeneracy is an approximate result as the [0, tt] field is 
small enough, (iv) In the [tt, 0] and [tt, 7r]-field cases, we 
analyzed the magnetizations and the magnon dispersions 
fFigs. [TlllSII . The inhomogeneous magnetization in the 
ladders were predicted. Moreover, it was shown that the 
[tt, 0] and [tt, tt] field do not induce any critical phenom- 
ena at least for the weak rung-coupling regime. This is 
in contrast with the the gap reduction by the [0, tt] field. 

The new even-odd nature and the quantum phase tran- 
sition between two critical phases in the [0, 0]-field case 
are most fascinating among all the results. However, one 
has to remember that our NLSM strategy is originally 
based on the case without external fields. Therefore, 
within such a strategy, one can not essentially provide 
a quantitative prediction for the case where the uniform 
field is so large that magnons are condensed. We will re- 
visit the magnon-condensed state in the frustrated tubes 
using other methods elsewhere. Furthermore, we will 
discuss half-integer-spin ladders and tubes in the near 
future. 

Besides our frustrated spin tube, (as we already stated) 
other mechanisms generating the magnetization cusp 
have been known However, such mecha- 



nisms usually require artificial or fine-tuned interactions 
in the models. On the other hand, the structure of spin 
tubes is quite simple, and it was shown in Sec. lIII B3l that 
the cusp in the tube is stable against some perturbations. 
Thus, wc think that our scenario of the cusp has a higher 
possibility of realization compared with other ones. 

Our previous work)i4 based on the perturbation theory 
and bosonization techniques, shows that the 2-leg spin-S* 
AF-rung ladder with the [tt, 0] field has 2S critical curves 
in the sufficiently strong rung-coupling regime, and they 
vanish in the weak rung-coupling one. This prediction is 
consistent with our analysis for the weak rung-coupling 
case. Both studies, however, cannot explain how 25 crit- 
ical curves fade away. 

It is worth noticing that all staggered ([0,7r], [7r,0], 
and [7r,7r]) fields generally make triply degenerate spin-1 
magnon bands split into the doubly degenerate trans- 
verse modes and single longitudinal one within our SPA 
framework. It has already known^itSSiSi^ that the same 
type of the band splitting appears in the spin-i AF chain 
with the staggered field, when the field is sufficiently 
small: the effective theory of such a spin-^ chain is a 
sine- Gordon model, which low-energy spectrum consists 
of the massive soliton, the antisoliton (these two are de- 
generate) and the breather (bound state of the soliton 
and antisoliton). Therefore, the band splitting of two 
and single ones may be a universal feature in ID AF spin 
systems with an alternating field around the isotropic 
[SU(2)] point. 
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APPENDIX A: SOME RESULTS OF SIMPLE 
MATRICES 

Here, we write down some results of simple eigenvalue 
problems, which are used in Sec. IIIII 

Let us define the following two N x N Hermitian ma- 
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trices appearing in Sec. Hill 
Ai Ai 

A ' 



V 



(Ala) 



B 



/Si B2 

B2 Bi 



V B2 



■■• A2 

A2 / 

B2 \ 



'■■ B2 
B2 / 



(Alb) 



Eigenvalues Am and corresponding eigenvectors dm of A 
are given by 



Am = Ai + 2A2 cos km , where km = 
— ^(^sinA;,„,sin2fc„i, . . . 



(A2a) 



N + 1' 

sin N km], (A2b) 



where m = 1 , . . . , and = 1 . Similarly, eigenvalues 
Bn and eigenvectors 6„ of B {N > 3) are 

2n7r 

2B2 cos kn, where k„ — ^ , 

2 



6« 



(A3a) 



'(^sin(fc„ + 7r/4),sin(2fc„ + 7r/4), 

. . . , sin(A^fc„ + 7r/4)) , (A3b) 



where n — q, . . . ,N — 1 + q [q E Z) and bf^ = 1. 
We next introduce a matrix 



Cii (^12 
C21 (^22 



(A4) 



where Cn is a normal matrix. The determinant of the 
matrix (|A4p satisfies the following well-known formula: 

det \C\ = det ICiil X det IC22 - C21C1VC12I. (A5) 



APPENDIX B: SINGLE CHAINS WITH THE 
STAGGERED FIELD 

We give a short review of the Green's function method 
for integer-spin chains with the staggered field, which was 
discussed in Ref. [s^ 

We start from the NLSM coupling a general external 
field J(x), which Euclidean action is 



S'e [n, a : J] = / dx 



Ce-SJ-' 



(Bl) 



where Ce is same as Eq. (|Kcjl . The effective theory 
for the staggered-field case corresponds to J — H/a = 



{0,0, H/a). Here, let us introduce a Green's function 
G°(x,x') as 



d^ + c'dl + 2igcX{^) G°(x,x') =(J'(x-x').(B2) 



After integrating out n, the action becomes 
5e[A: J]-^Tr[lnG"(x,y)] 

y dxdyJ(x).G"(x,y)J(y)+z j dxA(x). (B3) 

The SPE SSe[X ■ A/^Mj-.R^ed = ^ is evaluated as 
3G"(x,x) 

I dydzG"{y,K)G"{^,z)J{y) ■ J{z) - 1, (B4) 

which determines the saddle-point value Asp(x). One can 
represent several quantities using G'^ within the above 
SPA scheme. The staggered magnetization is 

< ^ 5(n"(x)) = « JSfs,:S\ 

^ ^ " (5J"(x) W"(x) 

- rfy[G"(x,y)+GO(y,x)]J"(y). (B5) 

The excitation structures are estimated from the singu- 
larities of real-time connected Green's functions. They 
are associated with imaginary-time (Matsubara) con- 
nected Green's functions through analytical continua- 
tion. The latter is 

= 52[(r,n(x)"n(x')^) - (n(x)")(n(x')'')] 
52 InZ 



^^^e[Asp : J] 

(5J"(x)(5J/3(x') ^ ,5J°(x)(5J/3(x') 



G"(x,x')+G°(x',x) 6^pI2 



/ dydz 



G°(x,z)G°(z,y) (B6) 



+G"(y,z)G°(z,x) ^iy){^ 



5jP{yL') J ' 



where the functional derivative 6^ /SASB means that first 
6/6A is performed, and then S/SB is done. The final 
term SXsp/SJ^ can be determined by the following trivial 
equation: 



= 



SSe 



(5J"(x') V5A(x) 
^A(y) 



. = dy 



■(5J"(x') ,5J"(x') V,5A(x) 



^A(x)a(y) 
SSe 



(B7) 



where S/SX\ j is the functional derivative under the con- 
dition that J is fixed, and S/SJ^lx in the final term 
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means the derivative with respect to the "expHcit" J"- 
dependence of SSb/SX\j. Through an easy calculation, 
Eq. l|B7p becomes 



(5J/3(x' 

G°(y, x)G"(x, x') + G°(x', x)G"(x, y 



(B8) 



where 

^(x,y) 



r(x,y) 



<5A(x)a(y) 
6r(x, y) + 25*2 J dzdwJ{z) ■ J(w) 

X G°(z,y)GO(y,x)GO(x,w) + 

G°(z,x)G"(x,y)G"(y,w)' 
GO(x,y)GO(y,x). 



(B9) 



Let us apply the above results to our staggered-field 
case, in which J = H/a. We assume that Agp is indepen- 
dent of X. It leads to the relation G"(x,x') = G'^(x — x'). 
The Fourier transformation of Go(x), therefore, can be 
defined as 



G"(k) 



dxe 



-ikx y^O 



G'^(x) 



^,(B10) 



where kx = kx — uinT. From Eqs. ljB4p and (jB10|) . the 
SPE is calculated as 

The final term denotes the deviation from the SPE 
of the no-field case. From Eq. ljB5|) . the staggered mag- 
netization is 



(B12) 



This result indicates that the staggered magnetization is 
parallel to H, namely = (0,0,7713). From Eqs. (|B6|I - 
(|B10|I . the Fourier components of the connected Green's 
functions are estimated as follows: 



G^^'2'2'(k) 



52G"(k) 



s'gc 



Gf (k) = G^^(k)-. 



3f(k) 



3r(k) + 2(777,/5)2G0(k)' 



(B13a) 
(B13b) 



where f(k) = ^^pGO(i[k + p])GO(i[k - p]) is the 
Fourier transformation of r(x). Other Green's func- 
tions all vanish. In order to know the excitation spec- 
trums, let us investigate the Fourier components of real- 
time Green's functions, G""(k) = G""{'k)\^^^-iz, where 
k = {z, k) and z = uj + irj. We refer the excitation 



modes determined from the poles of G^^'''^(k) = G^(k) 
[G^^(k) = G^(k)] to transverse [longitudinal] magnon 
modes. From Eq. (jB13a|) . the imaginary part of G^(k) 
is 



amGf(k) 



2e{k) 



S{uj ~ e{k)) ~ Situ + e{k)) . (B14) 



The delta-function singularity means that the transverse 
modes are exhausted by the singie-magnon excitations 
with the dispersion exik) — (-(k). The transverse gap 
is defined as = €^(0). The mode e-rik) is doubly 
degenerate correspondingly to x and y directions, in the 
present SPA scheme. 

The singularity structure of G^ is much more involved 
than that of Gj. Here, we show only its results. For 
convenience, we introduce several new quantities, 

GH = G,^(k)/(5V), (B15a) 

G^'^^\uj) =me{3m)G{Lu), (B15b) 

ri(2)(t^) = 3[<Hepm)f (k)]/(2c,c), (B15c) 

M = m,/S, e = e(fc), (B15d) 

where we omit the indication of the wave number fc, and 
f (k) = f (k)l^ At T = 0, using (5°, we calculate 

as 



3mr(k) = 
inef(k) = 
where A{k,p') 
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4:A{k,p') 
dz 



es(cj2-4e(fc/2)2)sgn(cj), (B16a) 



2z3mT{k,z)V 



(B16b) 



\p'e{p' + k) 



^ (p' + mp')\ w = 

f v/(w2 _ 4g(fc/2)2)/(t^2 _ ^2p) _ ;,/2], Os is the Heav- 
iside's step function, sgn(ijj) is the sign function, and V 
means the Cauchy principal part. To derive Eq. ljB16b|l . 
we used the Kramers-Kronig relation. Using the new 
symbols (jBlSII . Eq. ljB13bp is simplified as 



Giuj) 



T{co) 



r(cj)(e2 - {u; + ir] f) + M^ 



(B17) 



Thus, G^'2(cl') are written as 



a(w)2 -f /3(a;)2 



(B18a) 
(B18b) 



G^H 
G^u;) 

where a(w) = (e^ - lu'^ + jf)T'^{uj) + 2r]Lur'^{uj) + AP 
and f3{uj) = 2r]UjT\uj) ~ (e^ - uj^ + r]^)r'^iuj). The pole 
structure of G^ gives the longitudinal mode eL{k). For 
< 2e(fc/2), in which F^ = 0, we have 



G\u) 



2rjuj 



1)— +0 



[e2 -Lj^ + 772 - A/2/ri(a))]2 + 4772cj2 

(B19) 



7rsgn(tj)(5[/(w)], 
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where f{uj) — uP' ~ ^ ~ M"^ jY^{Lo). From this, one sees 
that under the condition eL(fc) < 2e(fc/2), the longitudi- 
nal mode satisfies /(eL(fc)) — 0, i.e., 

eL(fc)2 ^ eT(fc)' + MVr\eL(fc)). (B20) 

If the lowest excitation of the longitudinal mode is lo- 
cated in fc = 0, the longitudinal gap is defined by 
Ai = eL(0). Under the condition < 2At, one can 
easily perform the integral in r^(AL)|fc=o- As a result, 
its explicit form becomes 



where y = ^t- From Eqs. HB20p and HB21|) . we can 
arrive in the equation fixing A^^SS 



= 1 
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1 — - arctan 



(B22) 



Zg arctan 



2^A2,yy4^ 



at T = 0. On the other hand, does not have any 
(B21) singularities for |a;| > 2e(fe/2), in the SPA scheme. 
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